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The  purpose  o£  the  effort  reported  herein  was  to  develop  a 
mathematical  model  to  assist  in  the  design  of  laminated 
elastomeric  bearings.  The  approach  of  the  study  was  based 
upon  the  application  of  large-deformation  elasticity  and 
viscoelastic  mathematical  theories.  However,  the  complexity 
of  the  large- deformation  elasticity  theory  and  the  problems 
associated  with  computer  solutions  of  equations  of  linear  or 
classical  elasticity  for  the  particular  bearing  geometry  and 
for  materials  that  are  almost  lmcompressible  were  such  that 
this  analytical  work  was  discontinued. 

Although  the  effort  reported  herein  did  not  result  in 
attainment  of  the  program  objectives,  the  effort  did  result 
in  the  development  of  analytical  approximations  to  the 
application  of  linear  elastic  theory  and  yielded  results  that 
are  readily  applicable  and  very  close  to  exact  solutions 
for  bearings  for  light  loads. 

Future  efforts  by  this  command  toward  the  development  of 
design  formulae  for  laminated  elastomeric  bearings  will  be 
limited  to  an  experimental  investigation  approach. 
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ABSTRACT 


The  purpose  of  this  study  was  the  establishment  of  analytical  design  pro¬ 
cedures  for  laminated  elastomeric  bearings.  This  was  approached  with 
the  application  of  the  linear  mathematical  theory  of  elasticity  and  later 
with  nonlinear  large  Reformation  elasticity  theory. 

The  linear  theory  yielded  analytical  approximations  that  are  close  to  exact 
solutions  and  which  are  easily  applied  and  evaluated.  This  analysis  of  one 
typical  lamination  yields  the  distribution  of  stress  and  deformation  in  the 
elastomer  between  "rigid"  metal  lamina.  However,  the  limits  of  the 
linear  elasticity  theory  are  exceeded  for  greater  than  small  bearing  loads, 
indicating  the  need  for  the  application  of  the  more  comprehensive  large- 
deformation  elasticity  theory. 

The  large-deformation  theory  was  stated  and  the  equilibrium  equations 
were  derived,  but  the  solution  of  these  equations  was  not  carried  out. 
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PREFACE 


This  report  was  prepared  by  The  Franklin  Institute  Laboratories, 
Philadelphia,  Pennsylvania,  for  the  U.  S.  Army  Aviation  Materiel 
Laboratories  under  Contract  DA  44-  177- AMC-  1 10(T) . 

The  mathematical  studies  of  the  laminated  elastomeric  bearing  presented 
herein,  which  began  in  September  1963  and  were  concluded  in  December 
1964,  represent  the  sole  effort  of  The  Franklin  Institute,  for  which  Mr. 

R.  Clyde  Herrick,  Senior  Research  Engineer,  Applied  Mechanics  Labora¬ 
tory,  was  the  principal  investigator  and  author  of  this  report.  Acknowl¬ 
edgment  is  made  to  Mr.  T.  Y.  Chu  for  the  considerable  contribution  of 
the  analytical  approximations  and  computer  solutions  for  the  linear  theory 
of  elasticity  shown  in  Appendixes  A,  B,  and  C,  and  for  the  brief  presen¬ 
tation  of  pure  torsion  in  large  deformation  elasticity  shown  in  Appendix  D. 
Acknowledgment  is  also  made  to  Dr.  Barry  Wolf  for  the  many  formulations 
in  large  deformation  elasticity  leading  to  that  presented  in  Appendix  E  and 
for  his  contributions  to  the  work  in  linear  theory  done  in  conjunction  with 
Mr.  Chu.  The  contributions  of  Dr.  Kishor  D.  Doshi  and  Zenons  Zudans 
are  also  noted. 

This  constitutes  the  final  report  covering  the  first  application  of  the 
mathematical  theory  of  elasticity  to  the  problem  of  design  and  laminated 
elastomeric  bearings. 
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SUMMARY 


The  purpose  of  this  work  was  to  establish  a  design  procedure  for  lami¬ 
nated  elastomeric  bearings  based  upon  the  application  of  large  defor¬ 
mation  elasticity  and  viscoelastic  mathematical  theories.  However,  the 
complexity  of  the  large  deformation  elasticity  theory  and  certain  prob¬ 
lems  associated  with  computer  solutions  of  the  equations  of  linear  or 
classical  elasticity  for  the  particular  bearing  geometry  and  for  materials 
that  are  almost  incompressible  were  such  that  this  work  could  not  be 
completed. 

In  the  course  of  this  effort,  analytical  approximations  to  the  application 
of  linear  elasticity  theory  were  developed  that  not  only  yield  results 
that  are  very  close  to  exact  solutions  but  are  easily  applied  and  evalu¬ 
ated.  However,  these  approximations  are  for  linear  elasticity  theory, 
the  limits  of  which  are  exceeded  for  greater  than  very  light  bearing  loads. 
While  not  quantitatively  applicable  for  bearings  under  full  loads,  they  do 
give  much  information  concerning  elastomer  behavior  under  light  loads. 

The  equilibrium  equations  for  large  deformation  elasticity  theory  in 
cylindrical  coordinates  were  stated  for  combined  axial  compression  and 
torsion,  but  were  fully  derived  for  the  case  of  compression  only. 
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CONCLUSIONS 


While  the  objectives  of  the  study  were  not  reached,  much  useful  knowledge 
about  the  behavior  of  the  elastomer  in  each  lamination  under  load  was 
generated.  The  approximate  theories,  (a)  for  incompressible  materials 
and  (b)  for  almost  incompressible  materials,  are  useful  to  indicate  the 
relative  effect  upon  bearing  performance  of  bulk  modulus  (compres¬ 
sibility),  ratio  of  inside  diameter  to  outside  diameter,  and  width-to- 
thickness  ratio.  More  important,  much  was  learned  with  respect  to 
techniques  for  obtaining  a  useful  solution  of  the  nonlinear  large-deformation 
equations  that  is  necessary  to  predict  adequately  the  distribution  of  dis¬ 
placement  and  stress  within  the  bearing,  especially  along  the  bonded 
surface  under  realistic  loading  conditions.  The  original  beliefs  that 
large  deformations  are  present  within  the  elastomer  even  for  very  light 
compressive  loads  were  confirmed.  But,  although  much  is  known  as  to 
the  behavior  of  the  elastomer  under  load  from  the  application  of  linear 
elasticity  theory,  this  must  be  conditioned  with  the  statement  that  this  is 
qualitative  and  not  quantitative  for  full  loads. 
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INTRODUCTION 


During  the  experimental  development  of  the  laminated  elastomeric  bearing 
for  low -temperature  applications  (Technical  Documentary  Report  No. 
ASD-TDR-63-769,  prepared  for  Wright-Patterson  Air  Force  Base,  Ohio, 
November  1963),  the  need  for  a  mathematical  analysis  of  the  bearing 
which  would  form  the  basis  of  a  design  method  was  recognized.  Of 
primary  interest  were  the  magnitude  and  distribution  of  shear  and  normal 
stresses  between  the  elastomer  and  the  metal  lamina.  This  was  not  only 
to  enable  a  designer  to  specify  elastomers  with  adequate  strength  and 
adhesive  strength  (at  bond  to  metal  lamina),  but  to  determine  the  minimum 
thickness  of  metal  lamina  that  could  be  used.  Only  by  a  mathematical 
analysis  could  the  influence  of  compressive  load,  shear  load  (torsion), 
ratio  of  outside  diameter  to  inside  diameter,  width-to-thickness  ratio, 
and  shear  and  bulk  moduli  of  the  elastomer  be  studied. 

It  was  also  anticipated  that  the  deformations  and  strains  within  the 
elastomer  were  large  and  that  the  more  general  mathematical  theories 
of  elasticity  would  be  required  for  an  analysis,  especially  if  the  inter¬ 
action  of  compressive  load  and  bearing  torsion  observed  in  the  experi¬ 
mental  development  were  to  be  studied.  Linear  or  classical  theory  is 
made  linear  by  the  assumptions  that  both  deformation  and  strain  are 
infinitesimal,  thereby  removing  the  ability  to  study  the  interaction  of 
axial  compression  and  torsion. 

It  was  anticipated  that  in  this  study  the  state  of  stress  and  strain  in  the 
bearing  could  be  investigated  sufficiently  with  the  use  of  large-deformation 
elasticity  and  viscoelasticity  theories  and  that  a  computer  program  could 
be  devised  for  the  design  of  optimized  bearing  laminations  for  a  given 
application  and  loading  condition. 
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REVIEW  OF  EFFORT 


When  the  present  program  was  begun,  the  intention  was  to  go  directly 
to  large-deformation  elasticity  because  it  appeared  that  sufficient  theory 
was  at  hand.  The  intention  was  to  use  the  linear  or  classical  theory  only 
for  guidance  and  to  establish  certain  "known  points"  in  order  to  check 
the  large  Cicformation  work. 

Soon  after  the  large-deformation  formulation  was  begun,  it  was  discovered 
that  a  simplification  could  be  made  if  it  could  be  said  that  plane  surfaces 
in  the  elastomer  parallel  to  the  bonded  boundaries  remained  flat  planes 
after  deformation.  Some  preliminary  work  with  the  linear  theory  of 
elasticity  indicated  that  these  surfaces  did  not  remain  flat  planes  at  the 
free  edge  where  interest  in  stress  and  deformation  is  greater,  and  so 
the  formulation  of  the  large-deformation  problem  was  made  more  general. 
By  this  time  it  was  apparent  that  the  equations  for  large  deformations 
would  be  made  up  of  many  terms  and  that  solutions  would  be  time  con¬ 
suming. 

At  this  time  during  the  program,  it  was  believed  that  an  investigation 
should  be  made  of  the  effects  of  small  amounts  of  compressibility  of  the 
elastomer  upon  the  displacement  pattern  and  consequently  upon  stress 
magnitude  and  distribution.  This  was  needed  because  the  only  consti¬ 
tutive  relations  (stress- strain  relations)  that  were  immediately  available 
were  the  "Mooney  Relations"  (Reference  5)  that  were  formulated  for 
incompressible  rubber-like  materials.  No  relationships  were  at  hand  for 
compressible  materials,  because  it  seems  that  most  rubber-like  mate¬ 
rials  are  used  under  such  loading  conditions  that  the  hydrostatic  component 
of  total  stress  is  relatively  small.  Hence,  volume  compression  is  small 
and  may  be  neglected. 

An  investigation  into  the  effect  of  small  amounts  of  compressibility,  using 
classical  elasticity  theory,  was  then  begun.  First,  the  literature  was 
searched  for  solutions  to  similar  elasticity  problems  (the  problems  of 
compression  of  a  cylinder  with  ends  fixed  for  no  radial  displacement). 

One  similar  problem  (Reference  6)  was  discovered,  and  although  the 
boundary  condition  was  that  of  no  radial  displacement  at  the  outside  edge 
only,  instead  of  along  the  whole  bonded  surface,  it  was  believed  that  this 
solution  from  the  literature  would  be  useful.  However,  although  it  was 
an  analytical  solution,  the  effort  necessary  to  extract  information,  even 
by  a  computer,  for  various  conditions  proved  to  be  considerable.  It  was 
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much  beyond  the  effort  considered  necessary  to  make  an  independent 
solution  of  the  equilibrium  equations  of  elasticity  for  the  real  boundary 
conditions  (a  bonded  surface)  by  means  of  a  digital  computer.  Therefore, 
this  latter  approach  was  started. 

The  first  choice  was  to  solve  the  equations  of  a  compressible  material. 
Appendix  A  shows  the  formulation  of  the  equilibrium  equations  and  the 
complete  computer  program  (in  FORTRAN)  for  the  iteration  procedure 
to  yield  deformations.  In  the  theory  of  elasticity  for  homogeneous  and 
isotropic  materials,  two  constants  are  necessary  to  describe  the  mate¬ 
rial  behavior.  Of  the  sets  of  two  constants  that  can  be  used,  G,  the 
shear  modulus,  and  v,  Poisson's  ratio,  were  chosen  because  the  value 
of  Poisson's  ratio  reflects  the  relative  amount  of  compressibility: 
v  =  0.  5  for  incompressible  materials,  and  v  <  0.  5  for  compressible 
materials.  Since  the  equilibrium  equations  in  terms  of  displacements 
contain  the  coefficient  \H-2v,  which  approaches  infinity  as  v  approaches 
0.  5,  it  was  decided  for  the  first  solution  that  a  moderately  compressible 
material,  say,  v  =  0.  35,  should  be  used.  The  value  of  v  could  then  be 
increased  to  as  close  to  0.  5  as  possible  in  later  solutions,  after  check-out 
and  proof  of  the  iteration  procedure.  Except  for  some  initial  trouble 
in  programming  finite-difference  methods  for  the  geometry  of  a  very  wide  , 
but  thin,  annulus,  the  iteration  procedure  worked  well.  However,  on  sub¬ 
sequent  trials,  the  iteration  method  would  not  converge  to  a  re  xsonable 
solution  of  Poisson's  ratios  greater  than  \>  =  0.  4.  After  this  was  worked 
with  for  a  period  of  time,  effort  was  stopped  on  this  computer  solution 
for  the  compressible  elasticity  theory,  and  all  furtner  investigation  of 
compressible  materials  was  accomplished  with  the  use  of  an  approximate 
solution  described  below. 

The  next  effort  was  the  formulation  of  the  equations  for  incompressible 
elasticity  theory  and  the  computer  solution  thereof.  This  required  a 
completely  new  formulation  because  the  coefficient,  1/1-2P,  is  infinite 
for  v  =  0.  5.  Consequently,  the  new  formulation  introduced  a  pressure 
term  (Reference  1,  page  79,  problem  4)  in  the  displacement  relations 
instead  of  a  volume  compression  term.  Appendix  B  shows  the  derivation 
of  these  partial  differential  equations  and  the  FORTRAN  program  for  the 
solution  by  iteration. 

The  iteration  scheme  proved  to  be  erratic.  While  effort  was  being  made 
to  improve  the  convergence  of  the  method,  an  effort  was  made  concur¬ 
rently  toward  an  approximate  solution,  also  shown  in  Appendix  B.  It 
was  not  until  some  results  of  the  approximate  theory  were  used  as  the 
starting  values  for  iteration  that  convergence  seemed  probable.  By  this 
time  considerable  effort  had  been  expended  upon  this  solution,  because 
the  approach  was  contined  inasmuch  as  each  innovation  introduced  seemed 
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to  improve  the  results.  The  last  computation  made,  of  which  partial 
computer  output  is  included  in  Appendix  B,  appears  to  have  been  converg¬ 
ing  as  planned,  although  the  convergence  was  very  slow. 

In  the  meantime,  approximate  solutions,  both  analytical  and  computer, 
to  the  compressible  cases  of  linear  elasticity  were  investigated.  It  was 
recognized  that  the  exact  solution  of  the  partial  differential  equations  of 
equilibrium  was  not  of  major  interest  in  this  study  but  was  only  to  guide 
the  formulation  and  solution  of  the  large-deformation  case.  The  search 
for  approximate  solutions  was  fruitful.  An  analytical  solution  based  upon 
the  linear  theory  of  elasticity  was  derived,  not  only  for  the  incompressible 
case  but  for  the  compressible  case  as  well;  that  is,  flexible  enough  to 
accept  Poisson's  ratios  in  the  neighborhood  of  almost  negligible  compres¬ 
sibility,  v  =  0.  495  and  v  =  0.  499.  These  values  represent  the  order  of 
magnitude  of  Poisson's  ratios  for  elastomers  of  interest  in  the  laminated 
elastomeric  bearing.  These  approximations  for  incompressible  and  com¬ 
pressible  theory  are  shown  in  Appendixes  B  and  C,  respectively. 

Appendix  C,  for  the  approximate  compressible  theory,  includes  an 
attempt  made  to  solve  equations  (C-17a)  and  (C-17b)  by  computer,  using 
the  Rung-Kutta-Gill  method.  The  flow  chart  and  computer  program 
(FORTRAN)  are  shown  at  the  end  of  Appendix  C.  The  computer  solution 
did  not  work,  however,  and  the  reason  was  discovered  from  the  analytical 
solution:  the  solution  consisted  of  e+x  and  e"x,  with  x  assuming  very 
large  values.  Inasmuch  as  this  approximate  solution  for  compressible 
materials,  as  shown  by  the  analytical  solution  in  Appendix  C,  was  ob¬ 
tained  by  applying  the  variational  theorems  of  linear  elasticity,  it  jatis- 
fies  equilibrium  approximately,  and  it  satisfies  the  stress  boundary 
conditions  on  the  average,  that  is,  across  the  thickness. 

Appendix  D  shows  the  derivation  of  equations,  based  upon  large-deforma¬ 
tion  elastic  theory,  for  the  case  of  torsion  only.  This  was  done  mainly 
for  investigation  of  the  large-deformation  problem  of  combined  com¬ 
pression  and  torsion.  Note  that  it  is  infinitely  less  complex  than  the  case 
of  large-deformation  compression,  as  shown  in  Appendix  E. 

The  real  stumbling  block  in  this  study  is  the  set  of  large-deformation 
equations  for  axial  loading  (Appendix  E).  These  involve  constitutive 
relations  (stress-strain  relationships)  that  are,  as  yet,  not  known  for 
the  compressible  case  although  they  are  generally  known  for  the  incom¬ 
pressible  case.  The  real  problem  is  the  size  and  complexity  of  the  highly 
nonlinear  equations  and  the  consequent  lack  of  assurance  that  a  solution 
is  the  right  solution.  This  is  the  basis  of  the  desire  to  establish  solutions 
for  light  loads  using  classical  elasticity,  the  establishment  of  a  check 
point. 
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RESULTS  OF  LINEAR  THEORY 


While  the  computer  solutions  did  finally  yield  satisfactory  results  for  the 
incompressible  case  and  for  the  compressible  case  where  Poisson's  ratio 
was  0.  4  or  less,  these  displacements  or  deformations  were  not  trans¬ 
lated  into  stress.  By  the  time  that  it  could  have  been  done  with  reliable 
computer  solutions,  the  analytical  approximations  were  available.  These 
yield  a  solution  that  is  closer  to  the  exact  solution  of  the  equations  of 
linear  elasticity  than  the  application  of  linear  elasticity  theory  is  to  the 
real  bearing.  Perhaps  it  should  be  mentioned  here  that  the  linear  theory 
of  elasticity  is  believed  to  yield  a  very  close  approximation  to  elastomer 
stress  and  deformation  for  light  loads,  but  only  for  light  loads,  as  will 
be  shown  presently. 

SOLUTIONS  FOR  INCOMPRESSIBLE  ELASTOMERS 


Results  of  the  computer  solution  (Appendix  A)  for  the  compressible  case 
are  not  shown  here  because,  for  practical  elastomers,  the  ratio  of  bulk 
modulus  to  shear  modulus  is  so  very  large  that  the  behavior  is  more  like 
that  of  an  incompressible  material.  While  the  computer  solution  in 
Appendix  A  worked,  it  only  worked  for  materials  with  a  Poisson's  ratio 
up  to  0.  40.  This  is  too  low  for  elastomers  of  immediate  interest. 

Since  meaninful  results  were  obtained  from  the  computer  solution  for  an 
incompressible  material  and  since  a  useful  approximation  was  also 
derived,  the  results  of  the  work  shown  in  Appendix  B  are  summarized 
as  follows: 

1.  By  examining  successive  iterations,  it  was  found  that  the  iterative 
procedure  as  programmed  in  Appendix  B  converges.  The  con¬ 
vergence  is  much  improved  if  the  iterative  procedure  is  started 
with  the  results  of  the  approximate  solution,  but  it  is  still  slow. 
This  is  evidenced  by  the  fact  that  the  computer  output  shown  in 
Appendix  B  is  from  the  600th  iteration,  and  certain  inconsis¬ 
tencies  are  still  present  at  the  inner  and  outer  boundaries. 
However,  this  disturbance  is  seen  to  exist  at  the  free  boundaries 
and  propagates  only  about  one-half  a  thickness  (2h)  into  the 
bearing  from  either  boundary;  therefore,  this  was  accepted  as 
a  good  solution. 
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2.  The  influences  of  the  boundaries  r  =  rj  and  r  =  rc  sre  manifested 
only  near  those  boundaries.  Away  from  the  bound,  ies  by,  say, 

5  thicknesses,  the  approximate  solution  and  the  computer  coin¬ 
cide.  Thus,  the  approximate  solution  can  be  used  to  predict  all 
regions  of  the  bearing  except  near  the  boundaries  r  =  r^  and  r  =  rG. 
These  can  be  investigated  by  computer  solutions  using  much 
funer  gridworks. 

3.  Within  the  limits  of  the  linear  theory  of  elasticity,  which  has 
been  grossly  exceeded  in  the  computer  example,  the  analytical 
approximation,  as  shown  in  equations  (B-7)  through  (B-23),  is 
a  very  accurate  method. 

Figures  1  through  4  present  data  gained  from  the  analytica1  approximations. 
Figure  1  shows  for  incompressible  materials  the  relationship  between 
average  pressure,  width-to-thickness  ratio,  and  ratio  of  axial  displace¬ 
ment  to  elastomer  thickness.  Thus,  for 

80  psi  (shear  modulus) 

250  (width-to-thickness  ratio) 


0.  002  (ratio  of  axial  deformation  to 
elastomer  thickness), 

125G  =  10,000  psi. 

However,  from  Figure  2,  it  is  estimated  that  the  associated  maximum 
shear  strain  is  f  rz  =  9u/9z  =  1.63.  This  strain  is  far  beyond  the  limits 
of  the  linear  theory  of  elasticity  because  the  linear  theory  is  based  on  the 
assumption  that  Qu/3z  is  very  small,  as  compared  to  1.0. 

From  Figure  3,  a  measure  of  what  is  meant  by  "light  loads"  associated 
with  linear  elasticity  theory  can  be  obtained.  Although  9u/3z  =  0.  2  is 
still  much  beyond  the  values  acceptable  within  the  linear  theory  of  elas¬ 
ticity,  it  is  not  sufficiently  large  to  change  the  order  of  magnitude  of  the 
results.  Assume,  then,  that  3u/9z  =  0.  2  is  accepted.  From  Figure  3, 
for  a  width-to-thickness  ratio  of  250,  Pave  =  15.5  G.  Thus,  if  G  =  80  psi, 
as  for  the  silicone  rubber  used  in  previous  programs,  then  average  pr<  s- 
sure  is  1240  psi.  A  bearing  of  this  geometry  and  these  materials  was 
tested  to  an  average  pressure  in  excess  of  40,000  psi.  The  rubber  did 
not  extrude  or  come  unbonded,  but  the  steel  lamina  broke.  Thus,  it  is 


then 


G  = 


r  -r. 
o  1 


2h 


w. 


ave 
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known  that  elastomers  with  this  geometry  are  capable  of  extreme  pres¬ 
sures  and  extreme  shear  strains. 


For  strains  that  are  not  excessively  large,  say,  e  rz  =  0.  2,  a  relationship 
between  average  pressure  and  the  maximum  shear  stress  as  derived  from 
the  incompressible  elastic  approximation  can  be  shown. 


When 


W 

250  and  —  =  0.00025, 


then 


and 


= 

erz  dz 


0.2 


Shear  stress  is 


then 


and 


o  =  2G  e  =  0.4  G  ; 

rz  rz 

max 


rz 


JEM  = 


ave 


0.4 

15.5 


rz 


=  0.0258  P 


max 


ave* 


This  provides  some  measure  of  the  relationship  between  shear  stress  and 
average  pressure,  although  the  limits  of  linear  elasticity  have  already 
been  exceeded  at  this  strain  level. 


APPROXIMATE  SOLUTION  FOR  COMPRESSIBLE  ELASTOMERS 


The  approximation  shown  in  Appendix  C  has  also  yielded  convenient 
expressions  for  the  investigation  of  stress  and  displacement  within  the 
elastomer  of  a  lamination.  These  equations  are  summarized  as 
follows: 


for  g  ,  use  equation  (C-47); 

zaverage 

for  a  use  equation  (C-48); 

z 

max 


9 


for  ratio 

°z 

o  ' 

use  equation  (C-49); 

and 

for  ratio 

Zavg 

°rz 

max 

0 

use  equation  (C-50). 

z 

avg 

Figure  4  shows  a  comparison  of  data  for  compressible  and  incompressible 
theory.  The  relative  position  of  the  curves  indicates  the  influence  of 
various  amounts  of  compressibility,  although  the  exact  placement  of  the 
curves  is  not  necessarily  accurate  because  they  have  been  calculated 
from  linear  theory  where  the  strains  exceed  those  allowable. 
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6  s  SHEAR  MOOULUS 

UJo  _  RATIO  OF  AXIAL  DISPLACEMENT 
h  “  TO  ELASTOMER  THICKNESS 


2.00  (IO*3) 


1.75  (IO-3) 


1.50  (|0“3) 


1.25  (IO*3) 


1.00(10"*) 


0.75  (IO-3) 


0.50  (IO-3) 


0.25  (IO-3) 


1  1 


hickness  ratio) 


ersus  Axial  Compression 
Le  Elasticity  Theory). 


— 

>—  V-  0.495 

— 

^  _ j — v-  0  499 _ 

V  =  '/2  (INCOMPRESSIBLE) 


V-  ^(INCOMPRESSIBLE) 


V  •  0  495 


V  =  0499 


(width-to-thicknoss  ratio) 


NOTE:  FOR  v  =  •/g  ,  THE  CURVE  REPRESENTS  Pmax  /Pavg  AT 
THE  BOUNDARY  OF  ELASTOMER  AND  METAL  LAMINA 


Figure  4.  Comparative  Data  for  Compressible 
and  Incompressible  Theory. 
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APPENDIX  A 


COMPUTER  SOLUTION  FOR  A  COMPRESSIBLE  ELASTOMER 


This  is  the  formulation  of  the  pure  axial  compression  problem  of  one 
lamination,  following  linear  (or  classical)  elasticity  theory,  that  was 
prepared  for  computer  solution  using  finite-difference  techniques  and 
iteration  schemes.  The  geometry  is  that  of  the  flat  thrust  bearing. 

Consider  a  thin  disk  of  elastomer  compressed  between  two  rigid  plates 
and  assume  that  the  disk  is  composed  of  a  compressible  material  with 
Poisson's  ratio,  v  .  The  shape  and  the  dimension  of  the  disk  are  shown 
in  Figure  5. 


j 

♦ 

— r—+ 

^ _ r 

+ — n — 

^  '0 

Figure  5.  Sketch  of  Coordinates  for  One  Lamination. 


Since  the  disk  is  circular,  it  is  convenient  to  use  cylindrical  coordinates 
(r,  0,z).  The  symmetry  property  of  this  problem  indicates  that  stresses 
and  displacements  are  functions  of  r,  z  only. 

The  equations  of  equilibrium  are  thus  simplified  to 


and 


a?w  |  2(1  -»)  a2w  I  1  a2ji 

ar2  1-20  sz2  1-2“ 


I  22  + 

r  dr 


1  du  _ 
r  dz 


(A-  la) 


1  d2w  .  Sji  ,  2(l-v)  d2u  ,  2(l-u)  1 

l-2v  drdz  -  2  l-2v  -  2  l-2v  r 
dz  dr 


=° 


(A- lb) 
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where  u  and  w  are  the  displacements  in  radial  and  axial  directions, 
respectively. 


The  boundary  conditions  are 


U  =  0, 

u  =  0, 

°r  =  °> 


w  =  -6  at  z  =  +h 

w  =  6  at  z  =  -h 

a  =  0  at  r  =  r  and  r  =  r . 

rz  oi 


(A-2) 


Since  the  mid-plane  (z  =  0)  is  a  plane  of  symmetry,  the  boundary  con¬ 
ditions  can  be  replaced  on  z  =  -h  by  the  symmetry  conditions  at  z  =  0. 


At  z  =  o,  w  =  0 


and  ~  =  0  ; 
dz  ’ 


(A-3) 


thus,  the  problem  can  be  solved  only  for  the  domain  z  =  0  to  h  and  r  = 


ri  to  ro- 


Because  of  the  complexity  of  this  problem,  a  computer  program  is 
written  to  solve  this  set  of  simultaneous  equations.  The  equations  of 
equilibrium  and  the  boundary  conditions  are  written  in  finite-difference 
form  as  follows: 

wi-l, j  +  ^  +  2R(l)]  wi+l, j  +  21-2u^  wi, j-1 

+  wi , j+1  +  (*)(dS)(ui+l,j+l"  Vl,j+1  +  Vl,j-l“  ui+l,j-l^ 

+  (^iiy)(i^)(ui,j4-i  ~  wi,j =  0 


and 


*  ^l-2v^Wi+l,j+r  wi-l, j+1  +  wi-l, j-1  "  wi-l,j+l^  +  ui,j-l 


(-\  4-  - ^  ^  4.  2(l-v)  - b 


i,  j+1  l-2v  (1  2R(l)  ui+l,j  l-2u  (l  “  WTJ*  ui-l,j 

2 

M  J  I  I  —.11  I  .  - 

(A-  4b) 


-  2  [2=432  + 
*  l1-2u 


2R2(I) 


( 


4^)^.  •  =  o. 

1-2U  'Jui,j 


18 


For  a  general  point  (i,j)  in  the  body,  and  boundary  conditions 


U-  -  0 

i,mmm 

w-  =  -6 

i,raram 


Wi,2  =  0 

ui,3  =  ui,l 


z  =  h 


( A  -  5 ) 


>z  =  0, 


(u0  .  “  U .  )  +  7^77  (fr)u„  ,  +  7^7“  (Wo.^-W  n  )  o 


r  =  r 


r  =  r 


i=iirlu3,j'ui,jy  i-2p  V  2>j  i"2'’  2^fi  \ 

1  VU2,J+1“U2,J-1=-(W3,J‘W1,J) 

(  4^  (unnn+l,j  *  W,j>  +  Wfr  (f}  “nnn,^ 

+  l-2u  ^wnnn,J+l  "  wrmn,  j-1  ^  0 

[  “nnn,j+l  “nnn^j-l 


u - .x!  -  u —  ^  n  “  (wnnn+l,j  ”  wnnn-l,j^ 


MA-6) 


for  points  (i,  j)  on  the  boundary,  where  b  is  the  grid  size  and 

R(z)  =  ri  +  i  x  b. 

It  is  noted  that  since  the  material  is  compressible,  there  is  a  region  in 
the  body  where  the  material  experiences  only  hydrostatic  compression. 
This  knowledge  permits  further  reduction  in  the  domain  of  the  problem 
by  combining  the  two  end  conditions  into  or  e.  That  is  to  say,  the  set  of 
equations  is  solved  from  both  ends  (r  =  r^  and  r  =  rQ)  simultaneously, 
and  the  condition  of  hydrostatic  compression  is  used  as  a  boundary 
condition. 

Denote  ends  r  =  ri  and  r  =  rQ  by  k  =  1  and  k  =  2,  respectively,  and  define 
coefficients  as  shown  in  expressions  (A-7)  through  (A-9).  Then  simplify 
the  set  of  equations  to  those  that  appear. 
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Coefficients  for  the  computer  program  are 


A(l)  =  1-u 
A(2)  =  v/(l-u) 

A(3)  =  3-4U 
A (4)  =  l-2u 
A(5)  =  H/FM 
A (6)  =  l/(3-4v) 
A(7)  =  A(4)  x  A(6) 
A(8)  =  A(l)  x  A(6) 


k=l  R0(I,K)  =  RZER0  +  (FI-2.0)  A(5) 


C1(I,K) 

C2(I,K) 

C3(I,K) 

C4(I,K) 

C5(I,K) 

C6(I,K) 

D1(K) 


A(5)/R0(I,K) 

C1(I,K)  **  2 
1.0  +  0.5  *  C1(I,K) 
1.0  -  0.5  *  C1(I,K) 
A(3)  +  A(l)  *  C2(I,K) 
1.0/C5(I,K) 

+  1.0 


k=2 


R0(i jK) 
C1(I,K) 
C2(I,K) 
C3(I,K) 
C4(I,K) 
C5(I,K) 
C6(I,K) 
D1(K) 


R0UT  -  (RI  -  2.0)  *  A 
A(5)/R0(I,K) 

C1(I,K)  *  C1(I,K) 

1.0  -  0.5  *  C1(I,K) 
1.0  +  0.5  *  C1(I,K) 
A(3)  +  A(l)  *  C2(I,K) 
1.0/C5(I,K) 

-  1.0 


(A-7) 


(A-  8) 


(A-  9) 
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Note  that  in  the  following  expressions  the  double  sign  ±  or  =F  appears. 
The  upper  sign  refers  to  the  case  where  k  =  1,  while  the  lower  sign 
refers  to  that  where  k  =  2. 


Boundary  conditions  are 


WiJ  W3,J  ±  (U2,J+1  “  U2,J-1) 

ui,J  =  U3,J  -  (R(!i!k)  u2,J  +  W2,J+1~  W2,J-1 


(A-  1 0) 


Field  equations  are 


2(3-^) 


WI-1,J  +  ^  -  2^(1^)^  WI+1,J 


+  (2  +  £KwIjsKL  +  wj^ j+1)  ±  k  (1  +  ^^ui+i,j+i 

"  UI-1,J+1  +  UI-1,J-1  “  +  2R0(I,K)  (1  +  ft 


(A- 11a) 


04, ; ,  J  12 

4  *  2R02(I,K) 

"  WI-1,J+1  +  WI-1,J-1  “  WI+1,J-1^  +  UI,J-1  +  UI,J+1 


+  <2  +  ;> 


uI+1,J  +  (1  +  2R0(I,K))  UI-1,A 


(A-  lib) 


rv  aa-"  - 


Substituting  \/^  =  2u/l -2v, 


WI,J  =  2(3 -4v)  ^  *  2R*(I,K)^  wI-1,J  +  Cl  -  2R0(I,K)]  wI+1,J* 

+  >4v  +  wI,J+l^  -  2(3-4u)  ^uI+l,J+l“  uI-1,J+1 

-  b  (A- 12a) 

+  UI-1,J-1  "  UI+1,J-1^  +  2(3-4u)  ^2R^(l,K)^  (uI,J+1  "  UI,J-1^ 

ui,j  =  i  ^  “  T1  ^  ^  ^wi+i,j+i  “  wi-i,j+i 

C(3-4v)  +  (l-t>)  “j2 - 

+  WI-1,J-1  +  WI+1,J-1}  +  (1  -  2v)(uI,J-l  +  UI,J+1} 

+  2(l-v)  [(1  ±  2R0(I,k)^  UI+1,J  +  ^  ~  2R0(I,K)^  UI-1,J^*  ^  ^ 

From  the  definitions  of  the  coefficients, 

WI,J  =  °-5  *  A(7)  *  wi_i,j  +  °3(I.K)  wI+ljJ) 

+  A(8)*  (wj  „  +  Wj  J+1)  +  D1(K)*  0.125*  A(8)*  (ui+1)J+1 

-  '  ui+i,J-i'  +  0-25#  A(8)*  Cl(I,K)* 

(uI,J+l  ■  uI,J-r  (A- 13a) 

uI;J  =  0.5*  C6(I,K)*  D1(K)*  0.25*  (wJ+ljJ+1  -  V1)J+1 

-  »I+1,J-1>  +  A<A>*  <UI,J-1  +  \,m)  +  2* 

*(C3(I,K)  Uj+1  ,  +  C4(I,K)  Ul_1(J)).  (A- 13b) 

Boundary  conditions  ^end  equation)  are 

W1,J  =  W3,J  +  D1(K)*  (u2,J+l  ■  U2,J-1) 

u1> j  =  u3)J  +  D1(K)*  A(2)*  (2*  01(2, K)  u2)J  +  w2(J+1-  v»2> j^).  (A.  14) 
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Computation  of  stresses  is  as  follows: 


From  the  stress  -  strain  relation 


+  2Ge 


ij 


and  the  strain-displacement  relation,  stresses  can  be  calculated  at  any 
point  in  the  body  when  the  displacement  field  is  known.  Since  the  displace¬ 
ments  are  given  by  numerical  values,  finite-difference  methods  must  be 
used  to  calculate  strains. 


At  the  top  and  bottom  boundaries,  the  strains  are 


ee 


=  0 


€  =f* 

z  dz 


and  the  stresses  are 

0  =  (X  +  2G)  e 
z  z 


2G(1  -  v)  dw 
1  -  2v  dz 


a 

rz 


In  order  to  evaluate  the  derivatives  there,  backward  or  forward  differences 
must  be  used.  At  the  top  surface,  backward  difference  is  used. 


If  =  +  f(M*2)) 


Since  in  the  computer  result  u  and  w  are  nondimensionalized  with  respect* 
to  h,  A  z  =  h/m. 


For  m  =  4, 


^  =  2(3w  (I,mrarn)  -  4w(I,m)  +  w(I,m-l)) 


-  2(3u  (I,rarain)  -  4u(I,m)  +  u(I,ra-l)); 
oz 

az  =  "  4w(I,m)  +  w(I,ra-l) 

o  =  2G(3u  (I.nmm)  -  4u(I,m)  +  u(I,m-l). 
rz 
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COMPRESSIBLE  MATERIAL 

DIMENSION  A ( 8 ) *  K0(40l,2),  CI(4GI,2).  02(401, 2).  03(401,2), 

104(401,2),  05(401,2),  06(401,2),  01(2) 

COMMON  M,  KZERO,  ROUT,  H,  A,  RO,  N,  LITER,  ALPHA,  DELTA, 

1C0NST ,  K,  JITER,  I F I N,  Cl,  02,  03,  04,  03,  CA,  01 
110  READ  INPUT  TAPE  2,1,  KLUE  ,  ICLUE 

1  FORMAT  (215) 

GO  TO  (500, 100,400, 400), KLUE 

100  READ  INPUT  TAPE  2,  2,  N,  M,  R2ER0,  ROUT,  GNU,  H  .JITER, IFIN 

2  FORMAT  (214, 4E15. 8/216) 

READ  INPUT  TAPE  2,  3,  LITER,  ALPHA,  OELTA,  CONST,  PMU 

3  FORMAT  I  IS,  4E15.8  ) 

All)  ■  1.0  -  GNU 

A ( 2 )  -  GNU/AI 1 ) 

A ( 3 )  •  3.0  -  4.0-GNU 
A <  4 )  -  1.0  -  2.0-GNU 
FM  -  M 
A ( 5 )  >  H/FM 
A ( 6 )  -  1.0/AI3) 

AIT)  ■  A  (4 ) «A( 6) 

A ( a  >  -  A( 1 ) -A t  6) 

CALL  POUMP  ( A , A( 8 ) , 1 ) 

K  ■  1 

20  CONTINUE 

GO  TO  (30,40) ,K 
30  NNN  -  N-2 

00  32  I  -  2, NNN 
FI  ■  I 

RO(l.K)  •  R2ER0  ♦  (FI  -  2.0)  •  A ( 5 ) 

Oil  I , K  )  -  A( 51/ROI I «K) 

C2( I ,K)  -  Cl(I.K)  •  C1(I,K) 

0  3  ( I ,  K )  -  1.0  ♦  (0.5  •  CKI.KII 
04 ( I ,K )  -  1.0  -  (0.5  •  Oil I.K)  ) 

C5( I.K)  -  A ( 3 )  ♦  All) *C2 ( I , K ) 

06(1, K)  -  1. 0/05(1, K) 

32  CONTINUE 
Dim  -  l.o 
GO  TC  50 
40  NNN  -  N-2 

DO  42  1-2, NNN 
FI  ■  I 

RO(I.K)  -  ROUT  -  ( F I -2 .0) • A ( 5 ) 

Cl(I.K)  >  A ( 5 ) /RO (I.K) 

02(1, K)  •  CKI.K)  •  Cl  (  I  ,K  ) 

C3II.K)  -1.0-1  0.5  •  Cl(I.K)) 

04 (I.K)  -  1.0  ♦  (  0.5  •  Cl ( I , K ) ) 

05 ( I , K  )  -  A ( 3 )  ♦  All) *02 (I.K) 

C6II.K)  >  1. 0/05(1, K) 

42  CONTINUE 

01(K)  -  -1.0 
GO  TO  60 

50  WRITE  OUTPUT  TAPE  10,10,  R2ER0.RUUT 

10  FORMAT  (1H1  3X  8HRZER0  >  F10.5,  3X  7HR0UT  -  F10.5//) 
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WRITE  OUTPUT  TAPE  10*11.  GNU,  M,  H,  N 

11  FORMAT  ( 3X  6HGNU  ■  F10.5,  3X  4HM  «  15.  3X  4HH  -  F10.5,  3X  4HN  ■ 
115//) 

WRITE  OUTPUT  TAPE  10.12.  LITER, ALPHA, DELTA,  CONST, PMU 

12  FORMAT  (IX  8HL1TER  •  15,  3X  8HALPHA  ■  F10.5,  3X  8H0ELTA  ■  F10.5 
1,  3X  8HC0NST  •  F10.5,  3X  6HPMU  ■  F10.5////) 

60  CALL  SOLUT 

GO  TO  189, 108), K 

89  GO  TO  (108,90)  ,  1CLUE 

90  K  *  2 

GO  TO  20 
108  CONTINUE 
GO  TO  110 
400  CONTINUE 
GO  TO  110 
500  CALL  OUMP 

END  ( 1 , 1,0, 0,0, 1,1,  1,0, 1,0, 0,0, 0,0) 
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SUBROUTINE  SOLUT 

0 IHENS  ION  4(6),  R0(40l,2),  01(401,21,  02(401,2.1,  C9(401.2). 

104(401.2).  05(401. 2) .  06(401.2).  01(2).  UI401.1*),  MI401.1C), 
1UUU(10).  HWM(IO).  UU(10)»  WW(10) 

COMMON  M.  RZERO.  ROUT.  H.  A.  RO.  N.  LITER.  ALPHA,  DELTA, 

1C0NST .  K,  JITER,  IFIN,  Cl.  02.  03,  04,  05,  06,  01 
MM  ■  M  ♦  1 
MMM  ■  M  ♦  2 
NN  -  N* 1 
NNN  -  N*2 
FM  ■  M 

00  10  J  ■  2, MMM 
U(NNN.J)  -  0.0 
FJ  «  J 

W ( NNN , J  )  >  -(FJ-2.0)*DELTA/FM 
10  CONTINUE 

00  15  I  •  2,NN 
U(I.MMM)  -  0.0 
M( I, MMM)  -  -OELTA 
Mil, 2)  •  0.0 
15  CONTINUE 

00  20  I  ■  2.NN 
00  20  J  -  2, MM 
FJ  •  J 

U(I,J)  -  -01 ( K ) *CONST*DELT A* ( 1  •  0— I (FJ-2.0)*A(5)/H)**2) 
1*EXPF(01(K)*(R0(2,K ) — R  0 ( l.K))) 

M ( I , J  )  ■  -( FJ-2.0) *0ELTA/FM 
20  CONTINUE 
IJ  -  1 

JJ ITER  •  JITER 
00  110  ITER  -  1, LITER 
00  25  I  >  2, NNN 
U(I,1)  ■  U( 1,3) 

Ml  1,1)  ■  -M( 1 , 31 
25  CONTINUE 

00  30  J  •  2, MM 

U ( 1 , J )  -  U(3,J)*01(K)*A(2)*( 2.0*01  (2,K)*U(2,J)*MI2,J*1)-M(2,J-1)) 
W < 1 , J  )  ■  W ( 3 , J ) +0 1 ( K ) • ( J ( 2 , J* 1 ) -U ( 2 , J- 1 1 ) 

30  CONTINUE 

U(l,l)  -  U ( 1,3) 

M(l,l)  «  -M( 1,3) 

DO  80  I  ■  2,NN 
00  50  J  ■  2, MM 

UEVAL  -  0. 5«C6( l.K) *(01(10  *0.25*  (Ml  1*1, J*l)-M(I-l,J*l)*W( l-l.J-l) 
1-M( 1*1 . J— l ) ) *A( 4 ) • ( U(  I.J-l ) *U( I.J+l)  )*2 .0* A  ( 1)*(C3I I , K ) *U( 1*1, J) 

1 *04 (I,K)*U(I-1,J))) 

UU(J)  *  ALPHA  •  UEVAL  *  11.0-ALPHA)  •  UII.J) 

IF  ( J-2)  45,40,45 
40  MM( 2 )  -  Mil, 2) 

GO  TO  50 

45  MEVAL  -  0.5*A( 7)*(C4(  I ,Ki*M( I  —  1 • J >  *  C3( I ,K) *MI 1*1 , J ) ) 
1*AI8I«IM(I,J-1I*MI1,J*1II*  DllK)*0.l25*A(B)*(UII*l,  J*l) 
1-U(I-1.J*1)*U( I-l.J-1 )-U( I*1,J-1))*0.25*A(8)*01( l,K> 
1*(U(I,J*1)-U(I,J-1)) 

MM( J)  *  ALPHA  •  MEVAL  ♦  (1.0-ALPHA)  •  Ml  I , J > 

50  CONTINUE 
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IF  (1-21  60,60,70 
60  00  65  I  ■  2, MM 
UUUIJ)  ■  UUU) 

WWW ( J )  «  WW(J) 

65  CONTINUE 
GO  TC  80 
70  00  75  J*2 » MM 

U( I- 1 , J  >  ■  UUU(J) 

HI I-1,J)  >  WWW! J) 

75  CONTINUE 

CO  79  J  ■  2, MM 
UUUU)  »  UUU) 

HHH(J)  ■  WH(J) 

79  CONTINUE 

80  CONTINUE 

00  85  J  «  2, MM 
UINN.J)  ■  UUUU) 

H I NN, J )  -  WWW ( J ) 

85  CONTINUE 

IF ( ITER  -  I F IN  )  101,  105,  105 

101  IF  (ITEM  -  JJITER)  110,104,104 

104  IJ  *  IJ*1 
JJITER  *  I J • J I  TER 

105  INCR  ■  1 
JJ  -  1 

KK  *  MMM  -  6 

102  IF  IKKI106,  107,107 

106  JJJ*  MMM 
GO  TO  108 

107  JJJ  «  INCR  •  6 

108  WRITE  OUTPUT  TAPE  10,1,  ITER 

1  FORMAT  (IH1  17HI TERAT ION  NO.  «  15///) 

WRITE  OUTPUT  TAPE  10,3 

3  FORMAT  (  7HU  ARRAY//) 

WRITE  OUTPUT  TAPE  10 , 2 , IIUI I , J ) ,  J-JJ,  JJJ) , I-l.NNN) 
WRITE  OUTPUT  TAPE  10,4 

4  FORMAT  I//7HW  ARRAY//) 

WRITE  OUTPUT  TAPE  10,2, I (W( I , J )  ,  J- J J , J J J ) , I « 1 , NNN ) 

2  FORMAT  1 61 4X  E15.8)  ) 

IF  I  J J J-MMM I  109,103,109 

109  KK  «  KK-6 

INCR  ■  INCR  ♦  1 
JJ  ■  JJt6 
GO  TO  102 

103  CONTINUE 

110  CONTINUE 
RETURN 

END! 1,1, 0,0, 0,1, 1,1, 0,1, 0,0, 0,0,0) 
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APPENDIX  B 


ANALYTICAL  APPROXIMATION  AND  COMPUTER  SOLUTION 
FOR  AN  INCOMPRESSIBLE  ELASTOMER 


Governing  equations  are  again  written  in  cylindrical  coordinates,  as 
shown  in  Appendix  A.  The  equations  of  equilibrium  remain  the  same, 
but  the  stress-strain  relations  become 


a.  .  =  P  6.  .  +  2Ge.  . 
ij  ij  ij 


(B  -  1) 


where  an  additional  unknown  pressure,  P,  is  introduced  into  the  relation 
due  to  incompressibility.  Furthermore,  since  there  is  no  volume  change, 


C. .  =  0 
li 


(B-2) 


and  strains  are  deviatoric.  Taking  into  account  the  axial  symmetry 
property,  the  nonvanishing  strains  are 


-  s  11 

e0  r 


e  =  r— 
z  oz 


€  =  l  +  §2£ 

erz  *  dz  Br  • 


(B  -  3) 


Substituting  the  stress- strain  relation  and  the  strain-displacement 
relation  into  the  equilibrium  equation,  we  obtain  the  governing  equation 
in  terms  of  displacement. 


Lji  ,  1  _  1 

Sr2  r  3r  r2  Sz2  G  - 


(B-4a) 


1  d*V  _  1  dp 


2  +  r  5r  \  2 
?  oz 


G  dz  *  *  *  * 


(B-4b) 
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Since  there  are  three  unknowns,  u,w,  and  P,  three  equations  are  needed; 
the  third  equation  is  furnished  by  equation  (B-2), 


€ii 


=  0 


or 


du 

dr 


Equations 


(B-4a,  b,  and  c) 

u  =  0 
w  =  -w 

0 

w  =  w 
o 


coupled  with  boundary  conditions 

z  =  +  h 

z  =  +  h 

z  =  -  h 

r  =  r.  and  r 

1  o 

r  =  r.  and  r 
l  o 


(B-4c) 


(B  -  5 ) 


completely  define  the  problem.  However,  the  deformations  are  extremely 
complex,  owing  to  the  imposed  boundary  conditions,  and  the  mathematical 
analysis  becomes  intractable.  Thus,  the  use  of  approximate  methods  and 
computer  solutions  becomes  necessary. 


APPROXIMATE  SOLUTION 

Because  of  the  geometry  of  one  lamination,  [  h«  (rQ  -  rj)  ]  ,  the  boundary 
conditions  arr  =  0  and  °  rz  -  0  may  be  replaced  by  an  average  condition. 


Let 


(B-6) 


that  is  to  say,  the  system  of  forces  acting  on  the  boundary  is  replaced  by 
an  equivalent  system.  By  Saint  Venant's  principle,  this  replacement  will 
produce  negligible  influence  at  some  distance  away  from  the  boundary. 


% 


Consider  a  solution  of  the  form 


u  =  (ar  -  br”1)  (l  -  z%*)  ... 
2 


2  /,  2  > 


w  =  -2az 


1  - 


3h  2 1 


(B-7) 

(B  -  8) 


where  a  and  b  are  constants.  Substituting  equations  (B-7)  and  (B-8)  into 
the  governing  equations,  it  is  found  that  equation  (B-4c)  is  identically 
satisfied  and  that  equations  (B-4a)  and  (B-4b)  yield 

P/G  =  a  2|  u-r  +  2a(l  -  +  C,  ...  (B-9> 

h  h  br 

where  C ^  is  a  constant. 


This  solution  satisfies  boundary  conditions  on  z  =  ±  h,  identically, 
provided 

a  =  f  ^  (B -  10) 

4  h  • 

The  average  boundary  condition  on  r  -  rj  and  rQ  gives 


2_ 

2h 


/ 


-h 


h 


2h 


r 

•J 


+ 

or 


2h 


8u  , 

a7dz 


2h 


(B  - 1 1 ) 


(B-  12) 


Substituting  equations  (B-7),  (B-8),  and  (B-9)  into  (B-ll)  and  (B - 12),  the 
last  equation  is  satisfied  identically  and  equation  (B-ll)  yields 


k 

3 


(a  +  br"2)  + 


b  l^r 


=  0 


or 


^2“  "  ^2  ^  r  +  8  br”2  +  f  a  +  =  0  .  (B-13) 

n  n  3  2 
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Z* 


!,Vf, 


Hence , 


2 

r. 

1 


a 


2 

r 

a  -a- 
2 
Yi 


2b  ,  ,4-2  8 

+  +  J  a  +  C2  =  0 


%  ^  ro  +  I  brn"2  +  I  a  +  C,  =  0 

v,*  o  j  o  3  2 


(B-  1  3a) 


a  h 


b  = 


_o 

2 


2\ 
r . 

A 


2U  r  /r.  +  % 

o'  i  3 


r2 

1 


O  J 


(B-  14) 


r  r2 

a  “2 

n 


^  U  r.  +  f  a  +  ^  b  r72 
h  1  J  J  1 


(B  -  15) 


The  pressures  at  the  top  and  bottom  boundaries  (bonded  surface)  can  be 
calculated  from  equation  (B-9). 


(B- 16) 


where  a  is  substituted  from  (B-10)  and 
b  and  C2  from  above. 

The  total  axial  force  on  the  bearing  is  computed  by  integrating  equation 
(B-  16)  over  the  surface  of  z  =  h  (the  bonded  surface). 


P/G 


z  =  ±  h 


a  ^2  ’  ^2  U  r  +  C? 
n  hr  2 
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F/G 


(P/G)  r  dr  d6 


=  2tt 


r 

o 


J 


r. 

i 


-  r-t^t 

n 


r  +  C2r 


dr 


F/G 


TT 


ar  2b 

2  “  2 
2n  li 


-U,  r 


n 


(B-  17 ) 


Average  pressure  on  the  bonded  surface  is  the  total  force  divided  by  the 
area,  or 


P 

ave 

G 


hu  r  +  +  C 

n  n  * 


r 

t  o 


(B  - 1 8) 


To  find  the  maximum  pressure  on  the  bonded  surface  (z  =  h),  radius  is 
found  at  which  pressure  is  maximufn  by  differentiating  equations  (B-17) 
and  then  substituting  that  radius  into  equation  (B-16). 


and 


1  dP  =  ^ar  _  2b  1 


G  dr 


2  r 


=  0 


r 


(B-19) 


Now 


P 

G 


+  C, 


(B-20) 
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The  stress  distribution  at  the  inner  and  outer  radial  boundaries  is  now 
approximated  by 


ri  =  22u  +  p, 

G  ^  dr  +  G 


(B-21) 


and  by 


where  from  equation  (B-7) 


&  =  L  +  /-,  z_  ] 

>+r*  (1_h ?l’ 


grz  _  dq 
G  ”  dz  ’ 


(B-22) 


where  from  equation  (B-7) 


Consider  now  a  typical  example: 

=  0.7  5  inch 

rQ  =  1.125  inches 

h  =  0.  75  (10'3)  inch 

(width-to-thicknes s  ratio  of  250). 

If  the  ratio  of  axial  compressive  deformation  to  the  elastomer  thickness 
is  chosen  as 

f  ■  U0->, . 

then 

a  =  0.75  (10"3) 
b  ~  1.156  (10“3)  h 
C2~  ~  1-415  d0"3)  • 
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The  pressure  on  the  bonded  surfaces  can  be  computed  from  equation  (B-9) 


P/G 


z=h 


a 


U  r  +  C 

li 


2 


-  0.75(10“3)  % 
h 


-  2.312(10"3)  Ur  -  1.415(10’3) 


and 


r 

r/h 

U  r 

0.75" 

103 

-0.2877 

0.825" 

1.1  (103) 

-0.1924 

0.90" 

1.2  (103) 

-0.1054 

0.975" 

1.3  (103) 

-0.0253 

1.05" 

1.4  (103) 

+0.0488 

1.125" 

1.5  (103) 

+0.1178 

P 


0. 0629(H)3  )G 
0.0915(103 )G 
0.0892(103)G 
0.058(103)G 
0 


Average  pressure  is  computed  from  equation  (B -18): 


P 

ave 

G 


r 


L  ....  1 

2 

far2  2b  ,  ,  b  ,  „ 

i?  2j 

- 

2  2 
,  r  -  r. 

r 

z=h  o  i 

- 

-i 

r. 

l 


=  -  0.0629  (10~3) 


or,  in  magnitude, 

?ave  =  0.  0629  (10"3)G  (psi  when  G  is  in  units  of  psi). 

The  total  axial  force  applied  is  the  average  pressure  multiplied  by  the 
cross-sectional  area  of  2.  2  inches^. 

F  =  0.  1384  (10^)  G  (pounds  when  G  is  in  units  of  psi) 

Maximum  pressure  occurs  from  equation  (B-19)  at 

r  =  VVa  =  0.931  inch 
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and  is  of  magnitude,  from  equation  (B-20), 

P  =  -0.  0941  (103)  G  (psi)  . 
max  '  r  r 

p 

The  ratio  of  — max  is  \  5 
*ave 


Radial  stress  on  the  inner  radial  boundary  is 


a 

r 


G 


2 


P 

G 


z/h 

a 

r 

0 

1.812  (10“3)  G 

1/4 

1.427  (10“3)  G 

1/2 

0.485  (10"3)  G 

3/4 

-1.176  (10“3)  G 

1 

-3.5  (10“3)  G 

surface) 


Shear  stress  on  the  inner  radial  boundary  is,  from  equation  (B-22), 

a  - 
rz  _  ou 

G  -  dz 


or 


o  =  0.812  r  G  (psi). 

rz  n  r 

This  example  is  also  a  test  of  the  approximation.  From  the  above 
calculations  for  stress  on  the  inner  boundary,  it  can  be  seen  that  radial 
stress,  ar,  is  of  the  order  of  (10-3  G);  and  with  G  =  100,  for  example, 
the  radial  stress  is  still  of  the  order  0.  1  psi.  This  is  small  as  compared 
to  the  pressure  developed  within  the  bearing.  Shear  stress  on  the  bound¬ 
ary  varies  from  zero  to  0.  812  G.  This  is  not  small  as  compared  to  other 
shear  stresses;  but  since  a  shear  stress  on  this  boundary  does  not  in¬ 
fluence  greatly  the  pressure  within  the  bearing,  this  is  not  a  source  of 
serious  error.  The  approximation  is  considered  to  be  very  good  for  a 
wide,  thin  annulus. 
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COMPUTER  SOLUTION 


Due  to  symmetry  with  respect  to  z  =  0,  only  h  '.If  of  the  bearing  need  be 
considered;  if  z  is  positive,  the  symmetry  condition  then  states 


w  (r ,0)  =  0 


s\ 


du 

dz 


=  0 

2=0 


dp 

dz 


2=0 


0 


w(r,0+)  =  -  w(r,0~) 


(B-  23) 


To  solve  this  set  of  equations  (B-4a,  b,  and  c)  by  the  finite-difference 
method,  it  is  convenient  to  find  an  equation  involving  P  explicitly.  This 
can  be  achieved  by  differentiating  between  equations  (B-4a)  and  (B-4b) 
rather  than  by  using  equation  (B-4c).  Thus, 


1  dP  d2P 

r  9r  +  ,  2 
dz 


0  .... 


It  is  noted  that  equation  (B-24)  implies 


(B-24) 


rather  than 


du  u  dw 
dr  r  dz 


=  0 


dr 


+ 


(B-25) 


(B-26) 


It  is  thus  necessary  to  specify 


du  ,  U  .  cM 
fr  +  r  +  te 


0 


on  the  boundary. 
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Aj 


i .  .1+1 

i-l.-i  _ 

i+l.-i 

> 

E 

i 

I 

1 

t 

Written  infinite-difference  form,  equation?  (B-4a),  (B-4b),  and  (B-24) 
become 


2  2 
2  fi  +  2l  )  +  Jl. 

2  U  +  c2  '  +  R( I 


D 


ui,j  "  (ui+l'j  +  +  Wi)  (i+l,j  “ 


,  2 


D 


(ui,j+l  +  +  2H  (Pi+l,j  “  Pi-1,^ 


(B-27a) 


I  D2' 

2  1  +  =r 

\  C2< 


wi,j  "  (wi+l,j  +  Wi-l,j)  +  2R(T7  (wi+l,j  "  Wi-l,j) 

D2  D2 

+  ^2  (wi,j+l  +  +  2C,H  (Pi,^l  “ 


(B-27b) 


1  + 


D_ 

C2/ 


D 


P.  .  =  (P.  -  .  +  P.  .  )  +  OdTtT  (P’_l7  •  ”  P"  T  •) 
if  J  i*Plf  j  i“l»J  2R(X)  i+lf  J  i”ifj 


+  c2  (Pi,j+l  +  Pi,j-1) 


(B-27c) 


whe  re 


P  =  P/G  j 
H  =  h  . 


u 


Ui,j  H  • 


Wi,j  "  H  *  and 
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Treatment  of  equations  at  boundary: 


Since  the  boundary  conditions  involve  derivatives,  fictitious  points  are 
introduced  outside.  The  values  of  those  points  are  found  from  boundary 
conditions. 


At  r  =  r. 

l 


W.  .  =  W0  .  -  U0  .  .  ,  +  U0  .  - 

l*j  3,J  2,J+1  2,J-1 


P, 


2,:  ui,j  '  u3> j 


2  ,  D  ,  , 

U,  .  =  .  +  rT/'r, T  .  +  -  (w  .  .  .  -  -  W_  .  .  ) 

1,:  3,J  R(Z)  2,J  C  2,j+l 


(B-28) 


At  r  =  r 


WNNN,j  WNE,j  +  UNN,j+l  "  ^N,  j-1 
PNN,j  =  ^E,  j  “ 

_ R _  D  , 

UNNN, j  ^EJ  "  R(NN)  '  C  lWNN,j+l  "  WNN,j-l) 


(B-29) 


These  equations,  coupled  with  field  equations  (B-4a)  and  (B-4b),  are 

enough  to  solve  for  the  values  at  boundaries  rj,  rQ  and  the  fictitious 

points.  However,  it  is  not  desirable  to  indroduce  any  fictitious  points 

with  P  as  unknown.  Hence,  -H  is  written  in  terms  of  forward  or  back- 

dr 

ward  differences,  and  equation  (B-4a)  is  modified  accordingly. 


2E  =  1_  (r  P 

dr  2D 


i+l»J 


-  3  P 


-  p  ) 

i+2,y 


(forward  difference). 


(B- 30) 


<£  l_  , 
dr  2D  ^ 


P.  .  -  4  P-  t  .  +  P.  _  .)  (backward  difference).  (B-31) 
1>J  i-l fj  i-2, y 


At  z  =  h,  u  and  w  are  given  but  P  is  found  through  equation  (B-24).  In 
order  to  achieve  this,  it  is  necessary  to  introduce  a  fictitious  point  for 
P.  The  value  of  P  at  this  point  is  found  by  using  equations  (B-4c)  and 
(B-4b). 
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This  set  of  finite-difference  equations  is  solved  by  using  the  successive 
relaxation  method.  Values  of  u,  w,  and  P  at  one  point  are  calculated 
from  values  at  neighboring  points;  they  are  ueva,  weva,  and  Peva.  New 
values  of  this  point  are  obtained  from  ueva,  etc.  ,  through 

unew  =  »  ueva  +  (1  -  a)  „«*»  . 

These  new  values  can  then  be  substituted  for  the  old  values  immediately; 
then  one  can  proceed  to  the  next  point. 

Note  on  computer  data: 

Owing  to  the  limitation  of  computer  memory,  the  dimension  of  the  bear¬ 
ing  has  been  changed  to 


r^  =0.75  inch 
(ro  "  ri)/2h  =  50 

_  3 

h  =1x10  inches . 


The  approximate  solution  is  used  as  the  first  guess. 


INPUT  =  N,  M,  RZERfif,  H,  PMU 
OELTA,  R0UT,  CONST,  LITER, 
JITER.IFIN,  ALPHA,  BATA 


l 


Figure  8.  Main  Flow  Chart. 


Figure  9. 


Solution  Subroutine  Flow  Chart. 
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LAMINATED  BEARING 


INCOMPRESSIBLE  MATERIAL 

D l MENS  I CN  RO( 500) ,  6(8),  AR1I500),  AR2I500),  AR3(500),  AR4I50Q) 
COMMON  M*  RZERO,  H,  RO,  N,  LITER,  ALPHA,  DELTA,  BATA, ROUT,  B,C,0 
I»IFIN»JITER»  AR1  ,AR2, AR3, AR4  , CONST 
110  READ  INPUT  TAPE  2/1,  KLUE 

1  FORMAT  (15) 

GO  TO  (500, 100, 300, 400), KLUE 

100  RE AO  INPUT  TAPE  2,2,  N,  M,  RZERO,  H,  PMU,  DELTA, ROUT, CONST 

2  FORMAT  (2I5,4E15.8/2E15.8) 

READ  INPUT  TAPE  2<3,  LITER*  JITER,  IFIN,  ALPHA,  BATA 

3  FORMAT  ( 3 15,2615.8 ) 

WRITE  OUTPUT  TAPE  10,10,  RZERO,  H,  PMU,  DELTA, ROUT,CONST 

10  FORMAT  (1H1  3X  8HRZER0  ■  F10.5,  3X  4HH  ■  F10.5,3X  6HPNV  *  F10.5, 
13X  8HDELTA  *  F10.5/3X  7HR0UT  •  F10.2,3X  8HC0NST  «  F10.5//) 

WRITE  OUTPUT  TAPE  10,11,  N,  M,  LITER,  JITER,  IFIN,  BATA,  ALPHA 

11  FORMAT  (IX  4HN  «  15, 3X  4HM  «  15,  3X  8HLITER  -  15,  3X  8HJITER  •  15/ 
13X  7HIFIN  «!  15,  3X  7HBATA  «  F10.5,  3X  8HALPHA  •  F10.5////) 

FM  -  M 

NN  «  N+2 

C  «  H/FM 

D  *  C*B AT  A 

B 4 1 )  -  0/C 

B( 2 )  «  B ( 1  )**2 

B43)  -  1.0  ♦  B( 2) 

B(4 )  -  i.O/C 
B( 5 )  •  1.0/0 
B ( 6 )  «  C*B( 1 ) 

B( 7 )  -  B<4)*0(1) 

HIS)  -  1 • 0/B ( 3 ) 

CO  150  I  •  2, NN 
FI  -  I  -  2 
R0( I )  •  RZERO  ♦  F1*0 
ARKI)  -  1.0/’  I) 

AR2 ( I )  *  D*A<*(4) 

AR3II)  -  2. 0*8(3)  ♦  AR2(  I  )**2 
AR4( I )  »  1.0/AR3(») 

150  CONTINUE 

CALL  POUMP  ( RO , RO ( 1 0 )  , 1 ,B, B( 8 ) , 1, AR1, AR1 ( 10) , 1 , AR2 ,AR2 ( 10) , 
U/AR3,AR3(  10 )  «  1  ,AR4 ,  AR4(  10)  ,  1 ) 

CALL  SOLUT 
GO  TO  110 
500  CALL  EXIT 
300  CONTINUE 
400  CALL  DUMP 

END ( 1,1 ,0,0,0, 1,1 / 1,0, 1, 0,0, 0,0,0) 
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SUBROUTINE  SOLUT 

0IMENSI0N  R0( 500 ) «  U( 500 , 1 3 ) ,W ( 500 , 13 ) ,  P(500,13) 

1 ,  B  (  8 ) ,  AR1  (  500 ) »  AR2 ( 500 ) .AR31500) ,AR4(500) 

COMMON  M •  RZERO*  H,  RO,  N,  LITER*  ALPHA,  DELTA,  BATA, ROUT,  B,C 
1,IFIN,JITER,AR1 , AR2 , AR3 , AR4  , CONST 
MM  ■  M*2 
NN  -  N*2 
MMM  ■  M*3 
NNN  «  N+3 
NE  «  N*1 
ME  *  M*1 
RRU  *  1.0 
RRW  »  1.0 
RRP  «  1.0 
AA  *  0. 75*OELTA 

BB*  AA* I  ROUT**  2  -  RZERO**2)  /  1 2 . 0» LOGF ( ROUT /RZERO)  -(4. 0/3.0)* 
11(1.0  /ROUT ) *• 2  -  (1.0  /RZER0)**2 ) ) 

DO  10  I  «  2,NN 
DO  10  J*2 » MM 
FJ  »  J  -  2 

U( I , J  )  «  (AA  •  RO ( I  )  -  BB  •  AR1 ( I ) ) •  (1.0  -(FJ*C)**2) 

W(I,J»  «  -  2.0  *  AA  •  ( F J*C)  *  (1.0  -  (FJ*C)**2/3.0) 

P ( I , J  >  *  AA*  ( R 0 ( I ) **2  -  RZERO* *2 ) *2.0*BB*  LOGF ( RZER0*AR1  (  I ) ) 
1*2.0  *  AA*  (1.0  -  (FJ*C)**2)  -  8.0*AA/3.0-4.0*BB/(3.0*RZERO**2) 
10  CONTINUE 

P ( 2 , MM)  «  0.0 
P I NN»MM )  *  0.0 
DO  20  I  ■  2, NN 
U( I ,MM)  *  0.0 
W( I, MM)  *  -DELTA 
W  (  I , 2  )  *  0.0 

20  CONTINUE 

DO  30  I  ■  2,NN 
W(I,1)  «  -W( I ,3 ) 

U(I,1>  *  U ( I ,3 ) 

P( 1,1)  ■  P ( I , 3 ) 

30  CONTINUE 
I NCR  -  1 
JJ»1 
KK-MM-6 

1123  IF  (KK)  1130,1131,1131 

1130  JJJ-MM 

GO  TO  1151 

1131  JJJ-INCR  *  6 

1151  WRITE  OUTPUT  TAPE  10,3 

WRITE  OUTPUT  TAPE  1 0, 2 , ( ( U( I , J) •  J*JJ»JJJ)»I*1 , NNN) 

WRITE  OUTPUT  TAPE  10,4 

WFITE  OUTPUT  TAPE  10,2, (  ( W(  I ,  J) ,  J-JJ, JJJ) , I-l.NNN) 

WRITE  OUTPUT  TAPE  10,5 

WRITE  OUTPUT  TAPE  10,2, ( f  PC  I , J) , J»JJ, JJJ) , 1-1  ,NNN) 

IF  (JJJ-MM)  1155,1160,1155 
1155  KK-KK-6 

I NCR  ■  INCR  ♦  1 
JJ  ■  JJ+6 
GO  TO  1128 
1160  CONTINUE 
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IJ  ■  1 

JJITER  »  JITER 

DO  200  ITER  *  1, LITER 

RU  *  0.0 

RW  «  0.0 

RP  -  0.0 

DO  35  J  *  2* ME 

U(l.J)  *  U(  3 , J  )  ♦  2.0*AR2(2)*U(2,J)  ♦  BI1) * ( W 1 2 , J+l I -W( 2 , J- 1 ) ) 
UEVAL  »  AR4 ( 2 ) * ( ( U ( 3  *  J  )  *  U(1,J)I  ♦  0.5*AR2 C 2 ) • ( U(  3, J )-U ( 1 , J ) ) 

1  ♦  B(2)*(U(2,J+1)  ♦  U(2,J-m  ♦  0.5*0* {  4.0*P  I  3  »  J  )  -  3.0*P(2,J) 
i-P(4*J) ) ) 

PEVAL  *  B  (  5  )  •  (  Ul  1 ,  J )  -UI3.JD 
IF  (  ITER  -  JJITER)  34,29,29 
29  UU  *  ABSFIUEVAL  -  U(2|JI) 

PP  *  ABSF I PEVAL  -  P  ( 2 »  J  )  ) 

IF  (RU  -  UU)  31,32,32 

31  RU  *  UU 

32  IF  (RP  -  PP)  33,34,34 

33  RP  *  PP 

34  U(2,J)=  ALPHA  *  UEVAL  ♦  (1.0  -  ALPHA)  *  U(2,J) 

P ( 2 , J  ) *  ALPHA  *  PEVAL  ♦  (1.0  -  ALPHA)  *  P(2,J) 

I F ( J-2 )  45,40,45 

40  GO  TO  35 

45  W(1,J)  =W(3,J)*B(1)*(U(2,J*1)-U(2»J-1)) 

WEVAL  *  0.5*  8(8)  • ( ( W ( 3 , J ) *W ( 1 , J ) ) *0 . 5* AR2 ( 2 ) 

1*(W(3,  J)-W(1,J) )*B(2)  « (W(2,J*1 )+W(2,J-l ) )+0.5*B(6) 

.1*  (  P (  2 ,  J  +  l  )-P  ( 2  ,  J- 1 )  ) ) 

IF  (  ITER  -  JJITER)  30,36,36 

36  MW  *  ABSF ( WEVAL  -  M ( 2 , J  )  ) 

IF  (RW  -  WW)  37,38,38 

37  RW  *  WW 

33  W ( 2 » J  )  -  ALPHA  *  WEVAL  ♦  (1.0  -  ALPHA)  *  W(2,J) 

35  CONTINUE 

DO  100  I  *  3, NE 
DO  80  J  *  2, MM 
I F ( J  -  KM)  60,65,65 

60  UEVAL  «  AR4 ( I )  *  ( ( U( I ♦ l , J ) *U<  1-1 , J ) ) *0. 5* AR2 ( I > * ( U( I *1 , J ) 
l-U( I-l.Ji DB(2)*(U(  I,J*1)*U(  I,J-1>  >*0.5*D*(P( 1*1 ,J)-P(I-1, J) ) ) 

PEVAL  «  0.5*K(8)*((P(I*1,J)*P(I-1,J))*0.5*AR2(I)*(P(D1,J) 

1-P(  I - 1 , J) )*B(2) *(P(  I, J*1 ItPil  ,J-l ) ) ) 

IF  (  ITER  -  JJITER)  64,59,59 
59  UU  *  ABSFIUEVAL  -  U(I,JJ) 

PP  *  ABSF ( PEVAL  -  P ( I , J)  ) 

IF  (RU  -  UU)  61,62,62 

61  RU  *  UU 

62  IF  (RP  -  PP)  63,64,64 

63  RP  «  PP 

64  U(I,J)»  ALPHA  •  UEVAL  ♦  (1.0  -  ALPHA)  *  U(I,J) 

P(I,J)«  ALPHA  •  PEVAL  ♦  (1.0  -  ALPHA)  *  P(i,J> 

IF  (J-2)  75,70,75 

70  GO  TO  80 

75  WEVAL  «  0.5*B(8>*((W(Dl,J>*W(I-l,jm0.5*AR2(I  )  *  ( W  (  I  *  1 ,  J  ) 

1  -W  (  1-1,  J)  )*B(D*(W(  1,  J+1DWI  I«J-D  )*0.5*B(6)*(P(  I,J*1)-P(  I  ,J-l  )  )  ) 
IF  (  ITER  -  JJITER)  77,74,74 
74  WW  «  ABSF ( WEVAL  -  W ( I , J  )  ) 

IF  (RW  -  WW)  76,77,77 
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76  KW  *  MW 

77  W(I,J)=  ALPHA  *  WE V AL  ♦  (1.0  -  ALPHA)  *  WU.J) 

GO  TU  80 

65  P ( I  * MMM )  =P(  I  ,ME)*4.0*8(4)*(W(1,MM)-WI  I  ,  ME  )  ) 

PEVAL  =  0 . 5*  H  (  9  )  *  (  ( P ( I+l ,MM ) *P (1-1 , MM)  ) *0 . 5 • AR2 II) • ( P (  I  ♦  1 ,  MM  »  - 
1PI I-l.MM) )♦  8( 2 )* (P ( I , MMM)*P( I, ME) ) ) 

IF  !  ITER  -  J J I  TER )  79,81,81 
81  PP  =  ABSF  «  PEVAL  -  P ( I , J) > 

IF  (RP  -  PP)  78,79,79 

78  RP  *  PP 

79  P ( I , MM )  =  ALPHA  •  PEVAL  ♦  (1.0  -  ALPHA)  •  P(I,MM) 

90  CONTINUE 
100  CONTINUE 

00  95  J=2,ME 

UINNN.J)  =  U(NE,J)  -  2.0»AR2(NN)*U(NN,J)  -  B ( 1 ) * ( W ( NN, J*1 ) 

1-W ( NN , J - 1  )  ) 

UEVAL  =  AR4(NN)*l  IU(NNN,J)  ♦  U(Nt,J))  ♦  0 . 5  •  AR2 ( NN ) • ( U ( NNN , J  ) 
l-U(NE.J))  ♦  9  (  2 ) * ( U ( NN , J*l )  ♦  U(NN,J-1)»  ♦  0.5»D»(  3.0*P(NN,J)  - 
. 4.0*P (NE , J )  ♦  P (N, J  )  )  ) 

PEVAL  =  B ( 5 ) *  <  U ( NE , J  )  -  U(NNN,J)) 

IF  (  ITER  -  JJ!TER)  95,86,86 

86  UU  =  ABSF (UEVAL-UINN, J) ) 

IF  (RU  -  UU)  82,83,83 

92  RU  *  UU 

83  PP  =  ABSF(PEVAL-P(NN»J)  ) 

IF  (RP  -  PP)  84,85,85 

84  RP  *  PP 

85  U ( NN , J )  *  ALPHA  *  UEVAL  ♦  (1.0  -  ALPHA)  •  U(NN,J) 

P ( NN , J )  =  ALPHA  *  PEVAL  ♦  (1.0  -  ALPHA)  •  P(NN,J) 

IF  ( J-2 )  89,87,88 

87  GO  TU  95 

89  W ( NNN , J !  =  W(NE , J)-B ( 1 )* (U(NN, J+l )-U(NN, J-l >  ) 

WEVAL  =  0. 5*8(8 )•( ( W( NNN, J) *W ( NE , J ) )+0.5»AR2INN) 

1*(W(NNN,J )-W(NE , J ) ) ♦B(2 )  • ( W( NN, J* 1 ) *W ( NN,  J-l ) )»0.5»B(6) 

l*(P(NN,J*l  )— P  (Nf4,J  —  1 )  )  ) 
lF  (  ITER  -  JJITLR)  93,89,89 
89  WW  *  ABSF(WEVAL-W(NN, J)  ) 

IF  (RW  -  WW)  91 ,93,93 

91  RW  =  WW 

93  W(NN.J)  *  ALPHA  *  WEVAL  ♦  (1.0  -  ALPHA)  •  W(NN,J) 

95  CONTINUE 

00  110  I  *  2 , NN 
U (  1,1)  *  U( 1 ,3) 

W( 1,1)  «  —  W ( 1,3) 

PI  1,1)  *  PI  1,3) 

110  CONTINUE 

IF  (ITER-IFIN)  120,125,125 
120  IF  (  I  TER- J J ITER)  200,135,135 
135  I  J* I J  +  l 


J J 1  TER 

l  «  I J* J I  TER 

RUU 

S 

RU/RRU 

RWW 

M 

RW/RRW 

RPP 

a 

RP/RRP 

RRU 

a 

RU 

RRW 

a 

RW 

RRP 

a 

RP 
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125  I  NCR  -  1 
JJ  «  1 
KK  ■  MM  -  6 
128  IF  (KK)  130.131*131 

130  JJJ-  MM 
GO  TO  150 

131  JJJ  «  I  NCR  •  6 

150  WRITE  OUTPUT  TAPE  10,1,  ITER 

1  FORMAT  ( 1 HI  1 7H I  TER  AT  I  ON  NO.  -  15///) 

WRITE  OUTPUT  TAPE  10,3 

3  FORMAT  (  7HU  ARRAY//) 

WRITE  OUTPUT  TAPE  1 0, 2 , 1 ( U( I ,  J I ,  J« J J , JJ J ) ,  I -1 , NNN ) 

WRITE  OUTPUT  TAPE  10,4 

4  FORMAT  (//7HW  ARRAY//) 

WRITE  OUTPUT  TAPE  1 0, 2 , ( ( W ( I  ,  J ) ,  J> J J  ,  J J J ) , I > 1 , NNN ) 

WRITE  OUTPUT  TAPE  10,5 

5  FORMAT  (//7HP  ARRAY//) 

WRITE  OUTPUT  TAPE  1 0, 2 , ( < P( I , J ) , J» J J , J JJ ) , I ■ l , NNN ) 

2  FORMAT  (6(4X  E 15. 8 ) ) 

WRITE  OUTPUT  TAPE  10,6,  RU,  RUU,  RM,  RWW,  RP,  RPP 

6  FORMAT  (//2X  5HRU  «  E12.6,  2X  6HRUU  >  E12.6,  2X  5HRW  «  E12.6, 
12X  6HRWW  -  E12.6,  2X  5HRP  «  E12.6,  2X  6HRPP  >  E12.6) 

IF  (JJJ  -  MM)  155,160,155 
155  KK  *  KK-6 

I  NCR  -  INCR  ♦  1 
JJ  -  JJ+6 
GO  TO  128 
160  CONTINUE 
200  CONTINUE 
RETURN 

END( 1,1, 0,0, 0,1, 1,1, 0.1, 0,0, 0,0,0) 
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Note* 

1.  Upper  mh  letter*  on  grid,  Figure  S  refer  to  eoluans  on  so^uter  output. 

2.  Lover  o**e  letter*  refer  to  rov  froe  either  end  e*  noted. 

3.  U  err&y  denote*  red  lei  dlepleoeeent. 

•wple:  0.73755702  E  -  01  neene  U-  0.073755h 

4.  V  erreiy  denote*  ex  1*1  dlapleoeeent. 

ejou^le:  0.43044931  E  -  02  wen*  v  ■  0.004304h 

5.  P  errey  denote*  preiaur*. 

•nkjqtle:  0. 13768537  E  -  00  aeen*  P  -  0.137680 


Figure  10.  Geometry  and  Key  to  Computer  Output. 
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HI**!  ION  NO.  •  *00 


UliWi  «  MM1  , 

— v  Column  0' 


m —  71 1148112-01 

JJ®"  ■  ^-0.6107128)1-01 

nowe  *’^L-o.*nu«i*«-oi 

**wi^-o.657io)oii-oi 
-0.44)09)522-01 
-0. *4*150542-01 
-0.61)591142-01 
•0. 57M2  77*2 -01 
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APPENDIX  C 

APPROXIMATION  FOR  A  COMPRESSIBLE  ELASTOMER 


Consider  a  thin  disk  of  rubber  bonded  between  two  rigid  plates,  the 
rubber  being  assumed  as  compressible. 


The  dimensions  of  the  bearing  are  denoted  by 

ri  -  inner  radius 
r0  -  outer  radius 
2h  ~  thickness. 

Owing  to  the  difficulties  encountered  in  solving  the  equations  of  equilib¬ 
rium,  an  approximate  solution  through  variational  methods  will  be 
derived.  In  this  case,  the  theorem  of  minimum  potential  energy  is  used. 
A  displacement  field,  u,  w,  satisfying  the  displacement  boundary  con¬ 
dition,  is  assumed  and  put  into  the  potential  energy: 

U  =  i  Jy° ij  eij  dV  .  (c-l) 

The  theorem  of  minimum  potential  energy  states  that 

6U  =  0.  (C-2) 

Thus,  by  taking  the  first  variation  of  the  expression  U,  two  differential 
equations  and  the  associated  natural  boundary  conditions  are  obtained. 
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Assuming  that 


u  =  f(r)(l  -  ), 


w  =  ■jp  z  +  g(r)  z  (1  -  ^) 


|  =  f'(r)(l  -  £) 


I  -  «'W  » &  -  %> 

n 


£  -  *  +  «(r)(l  -  3 


m  =  £1x1  (i  _  a_). 


o  r 


Let 


4  -|a  +  »  +  |H  . 

dr  r  dz 


Then 


<L  4  e,  4  =  W  =  XA2  +  2G  [e  e  +  eQe0  +  e„e„  +  28  i] 


'ij  iJ 


r.  u  a  r,  a  w  w  AiW  w 

r  r  90  z  z  rz  rz- 


where 


and 


4  -  (f '(r)  +  f)(l  -  4  +  +  g(r)(l  -  3  4> 

h*  h 


o  #»  2  2  2  T.r~2  ~  2 

42  =  [f'(r)+f]  +T-  +  g2(r)(l  -  3 

h  h  h 


+  2[f  '(r)+  f]  (1  -  4^)  +  2  £  g(r)(l  -  3  4 

h  h 


+  2[f(r)+i]  g(r)  (1-4  (1-3  4- 


LI  *3 


a  AM 


) 

& 

<?■ 


The  following  are  expanded  for  substitution  into  equation  (C-l): 


e  e 
r  r 


[f'(r)]2  (1  -  4)2 


,2  2  2 
ee*e  2  ^  "  .2} 

r  n 


V,  ■  )2  +  g^(r) (l  -  3  42  +  2(^)  .Wd  -  3  4 


e  =  iC-  £(r)  +  g'(r)  z(l  -  *-)] 


2  2 


rz  rz 


^  f2(r)  +  g-2(r)  z2(l  -  42  "  *  4  f(r)  6’(r)(l-  ^4 
h4  n  h  n 


rbtx-4*-- fhcL-»4-4j- 

J-h  h*  J-h  h  h 


■Ca¬ 


li  h 


2£  +  i2!]h 

3  h2  5  h4  -b 


-  2[h  -  |h  +  |h]  =  i|  h 


lh  (1  '  3  h^>dZ  =  U  ’  h^)]"h  “  ° 

/h  (1  “  ^)dz  =Cz  *  3  ^]-h=  3  h 
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r  *2  2  r 

|  (1  -  3  h)  dz  -  (1  - 

J  -h  IT  J-h 


2  > 

6  +  9  \)dz 

h  h4 


=  z  -  2  *5  +  %  ^]h 
h  5  h4  -h 

=  2(h  -  2h  +  |)  -  |  h 


l 


(l  -  *5)(i 
-h  n 


3  ’SJds 


i 

i 


2  4 

(1  -  4  -S  +  3  *r)dz 

-h  h  h4 


=  k  2?.  +  1  zl]  h 

3h2  5?“h 

=  2(h  -  *  h  +  |  h)  =  fj  h 


J-h*2**  =  ^]-h  =  3  h3 


r  ,  _2  2 

J  *2(l-¥ 

-/-h  IT 


■f 


0  4  6 
(z2  -  z  ^7  +  ’2t)  dz 
'-h  h4  h4 


-/lJ_2*£  +  ljZ)h 
■  V  5  h2  7  h4'-h 


2(lL  „  2  h3  +  I  h3}  =iLh3 
2l3  5h  7n;  105 


rh 

L 


*3)4  "  £  *3  -  4  h2  1 


-h 


2(ih3-ih3)=^h3 


56 


n*'1 


t 


rA2  dz  =  h[r(r)  +  “l2  +  2  +  |  h  g2(r)  +  |  h(^)[f(r)  +  f] 


+  h  Cf»(r)  +  “]  g(r) 


Jere r  dz  =  h  [f'(r)]2 

Je6eQ  dz  =  if  h  ^ 

Jezez  dz  =  2h(^)  +  |  h  g2(r) 


e 

rz  rz 


dz  =  ^3h  +  105  gf2(r)  ”  15  h  8*(r)]  • 


Equation  (C-l)  is  equivalent  to 

Jwdv  =  {XA2  +  X  (erer  +  e6c0  +  +  2erzerz))dV 

=  J[XA2  +  2G(er€r  +  e0e0  +  +  Ze^JlrdrdBdz 

-  2nj\  tjf  h  [f  >(r)  +  +  2  +  f  h  g2(r)  +  f  h(&)[f  >{r)  +  £] 

+  i|  h  [f  *(r)  +  g(r))  rdr 

+  2n jx  {j|  h  [f ,2(r)  +  ^]  +  2h  (*-§-)  +  |  h  g2(r) 

+  ^  f2(r)  +  h3g'2(r)  -  ^  h  f(r)  g»(r))rdr ;  (C-7) 


but 


WdV  =  0 


(C-8) 


« 
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and 

J *{||  h  [r* (r)  +  £]  [6  f'(r)  +  ]  +  y  h  g(r)  6  g(r) 

+  f  C6f»(r)  +  f]  +  h  [6f»(r)  +  jftg(r) 

+  if  h  '(r)  +  6g(r))rdr 

+J 20  Clf  6^'(r)  +  6f]  +  ^*h  g(r)  6g(r) 

+  ^  f(r)  6f(r)  +  ^  Y?  g’(r)  6g'<r)  -  yy  h  f(r)6g'(r) 

-  ^  h  g’(r)  6f(r)}  rdr  =  0.  (C-8a) 

Expanding, 

J\  [|f  h  [f'(r)  +  f]  +  f  h  (*|)  +  if  h  g(r)}  6f'(r)rdr 

*/»  «  h  f’(r)}  6f’(r)rdr 

+/ x  h[r(r)  +  f](i)  +  f  h(f  )1  +  If  h  ^}6f(r)rdr 
+ J 2G{||  h  ^  +  yh  f(r)  "  15  h  S 1  (r)}  6f(r)  rdr 

+J 2G{y~y  g’(r)  -  yy  h  f(r))  6g'(r)rdr 

+J X  {yr  h  g(r)  +  yy  h  [fT(r)  +  ^])  6g(r)rdr 

+ y 2G{yr  h  g(r)3  6g(r)dr  =  0.  (C-8b) 
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< 


From  the  first  line  of  equation  (C-8b), 

J \  h  [f'(r)  +  +  |  h  (^)  +  25  h  g(r))  6f'(r)rdr 

*r[ff  h[f >(r)  +  +  |  h(Sf)  +  ^  h  g(r)}«f(r)]^ 

“  \f r^lf  hI^:Crf ' (r) ]  +  f'(r)]  +  j  h(^)  +  25  h  |j;(ig(r)])4f(r)rdr 
“  Xrtif  h[f '(r)  +  +  |j  (&)  +  i|  h  g(r)}  6f(r)]'° 

-  if  iff  h  [f  "(r)  +  f  f  '(r)]  +  |  h(&)  J  +  i|  h  [g'(r) 

+  ^l]j6f(rjrdr> 

a 

From  the  second  line  of  equation  (C-8b), 

J2&  h  f  *(r))  6f 1  (r)rdr 

=  *r[ff  h  f(r»  6f(r)]^°  -  »J ^  k  ^  It  f '(r)]6f(r)dr 

=  2Gr(ff  h  f'(r))  tf(r)£°  -  2G  (ff  h)J[t”(r)  +  ^]6f(r)rdr. 


From  the  fifth  line  of  equation  (C-8b), 

[aljffe  h3  g'(r)  -  ~  h  f(r)}  6  g'(r)  rdr 

=  2GrCjj|  h3  g'(r)  -  ^  h  f(r)}  6g(r)]r° 

"  XJr  I?  h?  rg'(r)  "  15  h  f(r)1  6g(r)rdr  ‘ 
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Grouping  terms  in  equation  (C-8b),  we  have  the  differential  equations  for 

f  and  g: 

X(22  h  [f"  +  if '  -  %]  +  hg'}  +  2G(^f  Cf"  +  *r  - 


3hf  +  15h«'1=° 


(C-9a) 


X  l15  h  8  +  if  h  [f '  +  rf]1  +  *  (5^ 

+  fj  h  ff'  +if])  =0. 


^  r_n  +  A  »t] 

105  Lg  r  g  J 


After  rearrangement,  the  equations  become 

fH-fifl  _  A—  f  _  -5 - 2 —  f  +  i 

1  -  1  ^  1  2  X+2G  1  * 


2h 


\+2G 


g'  =  0 


(C-9b) 


(C-lOa) 


g"  +  “  g’  -  ^2  8  ~  Cr'  +  “  f]  =  0 . 


?ii 


2h 


2  G 


(C-lOb) 


The  appropriate  boundary  conditions  are 


|  h2  g»  -  f  =  0 


(C-lla) 
(C-  lib) 


at  r  =  r.  and  rQ. 


For  the  sake  of  convenience,  these  equations  are  nondimens ionalized 
with  respect  to  h. 


Defining 


f  = 


£ 

h 

L 

h 


g 


-«<£>, 


(C- 12) 
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0 


(C-  13a) 


the  set  of  equations  becomes 

fit  +  i  7f  _  A_  7  _  it  iz.^2  7+1  -i-  '<  i  => 

1  -  1  -2  1  U-D  1  ;  1-v  8 

r  r 

g”  +  ^  S'  -  =  _  2  —i -  (f »  + 1  f)  n  Oi 

s  -  8  i-2v  8  2  l-2u  u  -  '  1 

r  r 

and  at  rQ  and  r^, 

+  i«  +  f(g»)J-0 


Is'  -?-0, 


where  prime  denotes  d/dr ,  i)  -  Poisson's  ratio  and  the  following  re 
lation  is  used: 


_  2Gu 

1- 


In  the  following,  the  bar  will  be  dropped  from  quantities  b,  q,  etc.  , 
the  understanding  that  they  are  nondimens ional. 


Define 


a 


.2  1=2k 

4  1-v 


b  -  i 


1-v 


c  _  211=21 

c  l-2v 


d 


2  l-2v  * 


The  differential  equations  can  be  written  as 


* 

f ii  +  1  f i  -  i-  f  -  af  +  bg*  -  0 
r  r2 

gii  +  i  gt  -  Cg  -  d(f f  +  ^f)  “  0  • 


(C- 13b) 

(C-  14a) 
(C- 14b) 

(C- 15) 

with 


(C- 16) 

(C  -  17  a) 
(C-  17b) 
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Consider 


f  =  A  L^kr) 
g  =  B  IQ(kr). 

Where  I's  are  modified  Bessel  functions  of  the  first  kind, 

Ij/  =  k  I0(kr)  -  1  Ii(kr) 

I,"  =  k2  I  (kr)  +  ^2  I.  (kr)  -  ;  I  (kr) 

r 

I0'  -k  Ix(kr) 

I  "  =  k2I  (kr)  -  z  I.(kr). 
o  o  r  j. 

Substituting  into  equations  (C-17a)  and  (C-17b), 

[(k2  -  a)A  +  kbB]  I±  =  0 

[(k2  -  c)B  -  kdA]  IQ  =  0 
or 

k4  -  (a  +  c  -  bd)k2  +  ac  =  0. 


(C- 18) 


(C- 19) 


(C-20a) 
(C-20b) 
(C -21) 


It  is  noted  that  this  equation  is  the  same  as  equation  (C-24);  hence,  the 
same  kj,  k £  will  satisfy  this  equation  and 


dk. 

F  =  -  - —  e 

.  2  * 
ki  -c 


(C-22) 


Now  the  general  solution  of  equation  (C-17)  can  be  written  as 
f  =  A  I^k-jr)  +  B  I^r)  +  E  K^r)  +  F  K^k,!-) 


(C-23) 


dk,  dk,  dk, 

«  -  A  Jo(klr)  +  ~Z~  B  Vk2r)  -  E  TT  K0(klr) 


.  2 
k^  -c 


k2-° 


k1-c 


dk 


,  2 
k2"c 


F  Ko(k2r). 


(C-24) 
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The  derivatives  are 

f  =  Aft-jl^k-jr)  -  ~  I^k^r)]  +  B[k2IQ(k2r)  -  J  ^(k^)]  -  E[k 

t 

+  7  K^r)]  -  F[k2KQ(k2r)  +  J  K^r)] 


dk  dk  2  dk2 

g*  =  — f-  A  Ifar)  +  -f-  B  Ix(k  r)  +  E  Kl(klr) 
k^  -c  k2«c  k^-c 


+ 


F 


^(kgr). 


Thus, 

(k2  -  a)A  +  kbB  =  0 
-kdA  +  (k2  -  c)B  =  0  . 

In  order  to  have  a  solution,  it  is  required  that 


or 


bk 


k2-c 


=  0 


k4  -  (a  +  c  -  bd)k2  +  ac  =  0 


k2  =  J[(a  +  c  -  bd)  +  “N^a  +  c  -  bd)2  -  4ac] 
k  =  +  {£(a  +  c  -  bd)  +  £"\/(a  +  c  -  bd)2  -  4ac)^  . 


Ko(klr) 

(C-25) 


(C-26) 


(C-27) 


(C-28) 


(C-29) 


(C-30) 
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Since  Ip(kr)  =  -Ip(-kr),  it  is  noted  that  there  are  only  two  independent 
solutions  corresponding  to 


k^  -  {£(a  +  c  -  bd)  +  £  *\/(a  +  c  -  bd)2  -  4ac}^ 

and 

k^  =  {J(a  +  c  -  bd)  -  J“\/(a  +  c  -  bd)2  -  4ac}^  . 


In  both  cases. 


k.d 


(C-31a) 

(C-31b) 


(C-32) 


Consider  next 

f  =  E  Kx(kr) 
g  =  F  Ko(kr). 

Where  K's  are  modified  Bessel  functions  of  the  second  kind, 

KJ  -  -k  Ko(kr)  -  ~  Kx(kr) 

K"  =  k2  Kx(kr)  +  ^  Kx(kr)  +  |  KQ(kr) 

r 

-  -k  K^kr) 

Ko  "  k2  •c0(kr)  +  r  Vkr)‘ 

Substituting  into  equations  (C-17), 

C(k2  -  a)E  -  bkF]  ^  =  0 
[(k2  -  c)F  +  dkE]  Kq  =  0  . 


(C-33) 


(C-34) 


(C-35) 
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<r 


< 


Thus , 


(k2  -  a)E  -  bkF  =  0 


dkE  +  (k2  -  c)F  =  0. 


(C  -  36) 


In  order  for  this  set  to  have  a  solution, 

-bk 
dk 


(k2-a) 


(k2-c )J 


=  0* 


(C-  37) 


Substituting  into  equations  (C  -11), 

2  '“‘l2  2 

(y  “2 -  “  1)  CA  +  E  +  (“ 


k£-c 


+  F  K^r.)]  =  0 


dk. 


k2-c 


1)[B 


^  cOc 

(f  -±-  -  1)[A  I^k^)  +  E  K^k^)]  +  (f  -  1)[B  VVo* 

V0 

+  FK1(k2ro)]=0 


k2-c 


and 


i  i  dk, 

A[klI0(^lri)  "  77  Il^klri^  +  T-v  ^77  Vklri^  +  ^  .2  Vklri^ 


k~-c 


+  BCk^Ck^)  -  ^  Ii(k2ri^  +  l-v  ^7.  Il^k2ri^  +  ^  .2  Vk2ri^ 


dk. 


k2"C 


i  l  dlt, 

-  E[k1Ko(k1ri)  +  —  K1(k1ri)  -  ^  (77  ^(k^)  -  J  2  K0(kiri^ 


k£-c 


,  i  i  dk_ 

-  F[k2Ko'k2ri)+  r.Kl(k2ri)- 

1  1  k2-c 


w. 


4  1-u  h 


(C-  38) 
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i  ,  dk, 

A[kx  I0(kxr0)  -  i-  (-  +  i 


k£-c 


-.  ,  dk_ 

+  B[k2Io(k2ro)  -  r  Il(k2ro)  +  &  <r  Wo>  +  i  ~t  Wo» 


V> 


l  i  ok. 

-  ECkiWo* +  r  Ki(kiro)  -  *  <r  Ki(kiro)  - 1 7T  Ko(kiro)] 


k2-c 


F[k2Ko(k2r0)  +  ^  Kl(k2ro) 
0 


(—  K,  (k_r  ) 
1-u  r  1  2  o 

o 


=  .  £  (— 2_) 

4  uV  h  * 


^  ,2  Ko^k2ro^ 

k2"C 


(C-  39) 


For  very  large  values  of  kr,  the  modified  Bessel  functions  have  the 
asymptotic  behavior 


kr 

1  (kr)  ~  -  ! - 

H  J  2rrkr 


Ip(kr)~-! 


(kr-*00) 


(C-40) 


Using  this  and  a  set  of  new  constants  defined  by 


k,r  -i  k0r  , 

A»  =  e  1  1  A  B'  =  e  2  ■  - - B 


'2nk, r 
1  o 


2nkr 
2  o 


(C-41 ) 


"klri  1  ”k2ri  1 

E*  =  e  1  1  ■  ■  ■  ■  E  F»  =  e  — Z— :  F, 


(2  k 
TT  K2  i 
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the  four  boundary  conditions  are  written  as 


o  dk 

A' 

1  k^-c 


-75 


Mr\-r  )  dk  2  fr  k  (r.-r  ) 

1  7  4c  ^  ri 


O 

+  E '  (-  — — 

7,2 

k^-c 


n  dk 

1  +  F*  <f  -f- 

k2-c 


-  1)  -  0 


(C-42a) 


0  dk,2  0  dk  2  p  dk,  ... 


2 

_  e 


+  F'  ( 


2^L 

?  4= 


-75 


K  -k2(ro-ri) 

e  =0 


(C  -42b) 


^  kl(ri-ro} 

a 


i  k-^-c  v  i 


+  B'[k„  -  1-  +-!!- 
2  r.  1-v 


kd_  lr~  k^(r--r  )  ,  , 

+  i  “5^“)]  Me210  -  E»[k,  +  —  -  r2-  (— 

*  .  2  v/  r.  1  r.  1-v  r. 

k2«c  v  l  ii 


dk^ 

k-^-c 


)] 


,  dkn  c  „  w„ 

-F>[k2+i 

1  1  K2  C 


(C-42c) 
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o  k^—c 


-  E«[k,  +  —  -  (i-  -  £  -i-)]  /-i  e 

1  r  1-v  'r  2  2  'v  r 


r.  -k,  (r  -r. ) 

i  l  o  i 


o  k~-c  u  ‘o 


o  k~-c 


-  F<  [k. 


l  1  dk~  /r7  -k-Cr  -r.)  c  w  (C-42d) 

ro  1'U  ro  k2-c  ^ro  4llVh 

For  large  values  of  r0-rj,  A1,  B1  in  equations  (C-42a  and  c)  and  E',  F1 
in  equations  (C-42b  and  d)  can  be  disregarded. 

Hence, 


E»  =  - 


2,2 
7  k2~C 
2 

g  <*1 
7  k^-c 


i  i  k,  -c  l  i  dk_-c 


=  1  (JL.)  2b 
4  klV  h 


(C-43) 


2  ,2 

n  k  -C 


A'  =  - 


2  *L 

7  k2-c 


,  ,  dk.  ,  ,  dk_ 

*'  Oci-r  +  i^(r  +  ^>’  +  B-  £  (j^i-t2-)] 


kj-c 


k2-c 


4  vlV  h 


(C-44) 
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In  terms  of  A's,  the  general  solution  can  be  written  as 


fr~  k,  (r-r  )  fr~  k  (r-r  )  fr~  -k,(p 

=  A'  \J—  e  1  °+B  'V-re2  +E'\/—  e 


+  F'/re 


-k2(r-r0) 


dk. 


=  — — a* 

6  2  v  r 


kj-c 


k,  (r-r  )  dk. 


■I  \  *■  — *  _  /  q  /  r  k  _  ( r— r  ) 

e  0  +  B  '\J-r  e  4 


.2 

k2-o 


dk,  fr~  -k,  (r-r.  ) 


-  E* 


,  2 
k^-c 


7? 


dk0  /r\T  -k2(r-ri) 


1'*  i  -  F'  J—  e 
e  ,  2  v  r 

k2-c 


When  rQ  -  rj  is  large,  the  solution  can  be  further  simplified: 

(1)  Region  near  rQ, 

yr~  k,  (r-r  )  fr~  k_(r-r  ) 

f  e  1  °  +  B'Jf  e  2  ° 


%  ki(r-ro> 


«  =A'  ~2~'JTe 

k^-c 


dk„  /r~ 

+  *'-iLJ-r 


dk0  /rT  k2^r'ro^ 

,  2  v  r  e 
k2-C 


(2)  Region  near  r^, 


f  =  E>  s/I 


-k-^r-m) 


*  Y* 


r,  -k  (r-r.) 

e 
r 


g 


■  -  ^  ys  .•*. 


-k0(r-r^) 


k£-c 


.2  v  r 
k2-c 


(3)  Away  from  ends, 

f  =  0 
g  =  o  . 


(C  -  45) 


(C-46) 


(C-47 ) 


(C-48) 
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Calculation  of  stresses  is 


at  points 


thus,  at 


Average 


a  =  (X  +  2G)  Is*  +  x(*  + 
z  Sz  Sr  r 

0  =  G(&  +  I2) 

rz  Sz  Sr 


on  bonded  surfaces  (z  =  ±  h) 


O  at  bonded  surface  is 
z 


avg 


2  2  J  z 
z^-h  ro  "r  ri 


I 


z=+h 


_  )  ( - - - ) 

^  1-2U  M  2  2' 


rdr 


r  -r . 
o  1 


Sw  , 
—  rdr 
Sz 


/2G(l-v) \  r /wo \  4 

^  l-2u  '  Uh  '  "  2  2 

r  -r . 
o  1 


grdr}, 


(C-50) 


(C  -  5 1 ) 


* 


(C-52) 
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W  % 


since 


dk 


g  =  [A  I0(kir)  -  E  Ko(kir)] 


k^-c 


dk. 


+  I0(k2r) 

k2-c 


F  Ko(k2r)] 


fo 


ri 


grdr  =  { — | —  [Ar  I^(k^r)  +  Er  K-^Ck^r) ] 


kx  -c 


dr 
k2"c 


f-  [B  I-^k  r)  +  F  K-^r)]]  [  0 
-c  J  r 


-  ro  CA  Wo*  +  E  Kl(klro)]  +  ~2  [fiIl(k2r0) 

k^  -c  k2  -c 

+  F  KjOtjT^]} 


-  r.  [— f-  [A  ^(k^)  +  E  K^r.)]  +  — |—  [B 
k,  -c  ’* 


k2‘C 


+  F  K1(k^ri)  ]) . 


For  very  large  values  of  kr, 


e _ u 

ro  [“2 


[A 


kx -c 


+k,  r 
1  o 

e _ 

\/2nk-  r 
1  o 


-k,  r 
1  o 


+  eT= 


■]  + 


.  2 , 

r  k2  JC 
o 


CB 


k,r 
1  o 


V*Vo 


+  F 


]} 


7  1 


I 


I 


,  klri  -klri  ,  k2ri 

ri  [~Z~  CA  "ft . .  +  E  6  '  1  +  -4-  [B  -  — 

k^  -c 


V^Vi  Vn  klri  k;2  "C  \/2TTk2ri 


-k2ri 


+  F 


vfv7 


,  /FT  -k,  (r  -r . ) 

ro  t— |—  a-  +E'7re 

k,  -c  o 


k22-c 


[B' 


/FT  -k(r  -r. ) 

+  F'Me  0  1  ]} 
\J  r 
v  o 


i  k  2-c  V  ri 


+  -f-  CB’ 

k2'C 


Fo  •kl(ro-ri) 

slT^ 


+  F']] 


or 


ri 


*rdr  =  ro  (-2~ A' +  ri 


k1-c 


kl  c 


B’)  -  r,(— f— E'  +  — | —  F>) 
kx  -c  kx  -c 


and 


max 


=  (ii  . 

1-2U  vh  * 


(C- 53) 


z=h 


hence, 


z  max  _ 


w 

(— ) 


(C-54) 


'*  avg  £  (^)  Vgrdr 

r  -r.  J 
o  l 
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It  is  also  noted  that  the  maximum  shear  stress  occurs  at  the  outer  edge; 
i.  e.  ,  r  =  rQ. 


Hence, 


0 

rz  max  _ 

a 

z  avg 


+  2G  f(rQ) 

2G(l-u)  A  L 

l-2v  'h  ”  2  2 

r  -r. 


_ f<r0> 

1=L_  A  — 4  _ 

l-2u  vh  “  2  2 

r  -r . 


(C-55) 


This  is  programmed  for  the  computer  as  follows: 

Let 

f  =  Xl  f »  =  ZL 
g  =  X2  g»  =  Z2. 

The  differential  equations  can  be  written  as 

(ZU'+izi-ij  +  ^  Z2  =  0 

r 

(XL) '  =  Z1 

(Z2),  +  i  Z2  -  X2  -  (  Z1  +J  XX)  =  0 

(X2) «  =  Z2  . 


Define 


DRAD  =  (RADN  -  RAD0)/FNR 
RAD  =  RAD0  +  FI#RAD  FI  =  I  -  1 
NRR  =  (NR/NP)  f-  1.0 

ND  -  number  of  intervals  skipped  during  printing 

(Zl)  * 

(XL) 1 

(Zl) » 

(X2)» 

Define 

DERZ1F  (A,B,C,D)  =  DRAD*|-|  +  +  (l.25*AA)#B  -  0.25*BB#C) 

DERX1F  (A)  =  DRAD*A 

DERZ2F  (A,B,C,D,E)  =  DRAD*  |  ^  #D+3.5*CC#(A  +  |)  j 

DERX2F  (A)  =  DRAD*A  . 

The  four  independent  solutions  of  XI,  Zl,  X2,  Z2  are 


XI 

S2(1,IJ), 

S2(2,IJ), 

S2(3,IJ), 

S2(A,IJ) 

Zl 

a(i,u). 

S1(2,XJ), 

S1(3,IJ), 

S1(4,IJ) 

X2 

S4(1,U), 

S4(2,U), 

S4(3,IJ), 

S4(4,IJ) 

Z2 

S3(l,IJ), 

S3(2,IJ), 

S3(3,IJ), 

S3(4,IJ). 

Where  the  first  index  denotes  number  of  independent  solutions, 

F(IJ)  =  A(l)  S2(1,IJ)  +  A(2)  S2(2,IJ)  +  A(3)  S2(3,IJ)  +  A(fc)  S2(4,I)  -  XI 
FF(IJ)  =  A(l)  S1(1,IJ)  +  A(2)  Sl(2,IJ)  +  A(3)  Sl(3,IJ)  +  A(4)  S1(4,IJ)  =  Zl 
0(IJ)  =  A(l)  S4(1,IJ)  +  A(2)  S4(2,IJ)  +  A(3)  SA(3,IJ)  +  A(4)  S4(4,IJ)  =  X2 
GG(IJ)  =  A(l)  S3(1,IJ)  +  A(2)  S3(2,IJ)  ^  A(3)  S3(3,IJ)  +  A(4)  S3(4,IJ)  =  Z2. 


=  Zl  +  -ir  XL  +  (1.25*AA)*(XL)  -  (0.25*BB)*z  2 
^  RAD2 

=  Zl 

-  1  afe  »  +  m  *Z2  +  3.5«cc*(zl  +  gg) 

=  Z2  . 
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Since  the  boundary  condition  states 


using 


and 


f '  +  (j^)  (f  +  i  I) 


w 


ri 


i  (_H_)  ^ 

4  vl-u'  h 


?'  +  (£>  &  +  i  i) 


w 


£  (^L)  -2 
4  U+t>'  h 


2/7  1 1  -  f  =  0 


2/7  g*  -  f  =  0, 


DD  =  GNU/ (l.O  -  GNU) 

W0  =  Wq/11 


F,  FF,  G,  GG, 


the  following  are  obtained: 

(si(ljl)  +  DD*  (s^lj'  +  0.5*S4(1,1))) 

+  (si(2,l)  +  DD*  +  0.5*S4(2,l)))  A2 

+  (si(3,l)  +  dd*(Smd/^  +  °-5*S4(3,l)))  a3 
+  (si(4,l)  +  dd*(~rad/^  +  0.5*S4(4,1)))  At  =  "  DD*w^*1.25 

(si(l,NRR)  +  DD*  +  0.5*  S4(1,NRR)))  Ax 

+  (S1(2,NRR)  +  DP*(S2^yRi  +  0.5*  S4(2,NRR)))  A2 

*(S1(3,HRR)  +DD*fi&gBl  +  0.5*  S4(3,NRR)))  A3 

+  (S1(4,NRR)  +  DD*  (^DN^  +  °’5*  s^>mR)))  \  =  "  DD*  1-25* 
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ft 


(|  S3(l,l)  -  S2(1,1))ai  +  (|  S3(2,l)  "  S2(Z,1))  A2 

+  (|  S3(3,l)  -  S2(3,1))a3  +  (f  S3(4,l)  -  S2(4,l))  A^  =  0 

(V7  S3(1,NRR)  -  S2(l,NRR)jA1  +  (2/7  S3(2,NRR)  -  S2(2,MR)^A2 
+  (2/7  S3(3,NRR)  -  S2(3,NRR))a3  +(.2/7  S3(4,NRR)  -  S2(4,NRR))A4 


=  0. 


Let  these  equations  be  represented  as 


{TRX}  [A]  =  [B]  . 


Solving  for  A, 


{TRX]-1  {Bj  =  [Aj. 
Since  B  can  be  written  as 


-  DD*  1.25*  w0  |0[  , 

ol 

{A}  =  -  DD*  1.25*  [TRX] 


-i/l 

SI 


Calculation  of  stresses  is 


0u  =  W6ij  +  2C  eij 


0  =  Xa  +  2G  e 
z  z 


-  (A+2G)e  +  X(cft  +  e  ) 
z  or 

=  U+2G)  f*  +  +2i). 

oz  dr  r 


a  =  2Ge 
rz  rz 


—  p  ( 1  5w\ 
_  G  +  3^- 
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Consider  points  on  flat  surface  (z  =  ±  h), 


&  =  a  =  o 

dr  r 

|n-0- 

dr 

Thus, 

°z  =  |f  =  (^+2G)  [j^2  -  2g(r)  ] 

=  (X+2G)  (w(Zf  -  2G)  =  2&=!!l  (wjf  - 

°rz  =  G  If  =  -  20  ^  “  *  20  *  • 
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_  e:. _ 


2G) 


sarr  -•  -*•**-* 


INPUT  =  GNU,  RAD0,  RADN, 
NR,  NR,  PMU,  W0,  H 


iffx 

i  r 

r-~ 


Fig  ure  12.  Subroutine  Int. 
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Fig  ure  14.  Subroutine  Stress. 
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AUTOMATH  SYSTEM  04/21/M  PAGE001 

•  •  •  SOURCE  PROGRAM  LISTING  ••• 

TITLCIEARING 

c  main  program  laminated  scaring 

COMMON  GNU.  RAOOt  RAON,  NR.  NP.  PMU.  MO,  H.X1 ,21 .) 2.22 . 

1S1.  SI,  SI.  S4,  P,  PP,  G,  GG,  SN,  ST 

DIMENSION  XI(4),Z1(4),22(4),X2(4),  S1(*«S1  )  *52  <*.91  I,SSI4.S1). 

154(4,51  ).  P (51) •  FPI51),  6(51),  GGISl).  !>N(51),  ST (51 ) 


10  RE  AO  4,GNU.RAOO,RAON,NR,NP,PMU,WO,H  1 

4  FORMAT  (SE15.S,  2I5/SE1S.I)  2 

IP  ENO  OP  PILE  200,201  S 

200  STOP  4 

201  PRINT  4,  GNU. R ADO, RAON, NR  5 

A  FORMAT  (1H1/SX.AH6NU  ■  P10.5,  SX.ThRAOO  •  P10.S,  SX.ThRAON  ■  P10.3  4 

1,  SX.5HNR  ■  15//) 

PRINT  T,NP,  PMU,  MO,  N  T 

T  FORMAT  ISX.5HNP  ■  If.  IX.4MPMU  ■  P10.5,  JX.SMPO  ■  P10.5,  • 

1SX.4HH  ■  P10.I////) 

20  CALL  I NT  5 

40  CALL  COE  10 

PRINT  )  11 

1  FORMAT  H1/4SX.20HD1SPLACEMCNTS  AND  DERIVATIVES//!  12 

PRINT  2  IS 

2  FORMAT  (//  14X.2HIJ  4X.11HP  FUNCTION  TX.llMG  FUNCTION  ISX,  14 

11SH  F  DERIVATIVE  TX.1SH  G  DERIVATIVE//) 

NRR  ■INR/NPI*!  15 

PRINT  S.UJ,  F(1J).  G(  IJI ,  FF(IJ),  GG(  I  Jl  ,1  J«1,NRR)  14 

5  FORMAT  (14X,|2,2X,ClS.a,SX.ElS.a,  14X,ElS,a,SX,ClS.ai  IT 

40  CALL  STRESS  14 

PRINT  101  15 

101  FORMAT  ( 1H1/4TX  2SH  STRESS  AT  FLAT  SURFACE)  20 

PRINT  102  21 

102  FORMAT  (//  ISX  2MIJ  Tx  14H  NORMAL  STRESS  10X  14H  SMEAR  STRESS)  22 

PRINT  10S,  ( I J.  SN(IJ),  ST(IJ) ,  IJ*1,NRR)  21 

10S  FORMAT  (SIX  12,  4X  C15.B,  5X  115.4)  24 

GO  TO  10  25 

END  24 
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AUTOMATH  SYSTEM  P01L301S  PAWOOl 


'nr' 


Ak  21?  ■  DER?  IF (21 t INT)  •  0.}*AKZlltXl(|NTl  •  0.5«AKXll.  22  <  INT ) 
WO.SaAKZZU  RAD*0.»«0RAtM 


•  »  o 

A  «t  « 


e 


DERXIFIAt  •  A  •  ORAD 

DERZ2F (A.BtC.D.E)  •  1-C/E  •  Zl.0»0/AA  •  ).9*CC*(A  •  B/E))»DrAD 


AUTOMATH  SYSTEM 


04/22/44 


PAGE003 


•  ••  SOURCE  PROGRAM  LISTING  ••• 

SUBROUTINE  1NV(A) 

C  ROUTINE  TO  INER  A  4  BY  4  MATRIX*  AC4»  4) 

DIMENSION  A(4t  4) 

00  9  I  •  1*  4 
B  •  All,  1) 

DO  1  K  *  1*  4 

1  A ( I •  K)  ■  A (I *  K)/B 
DO  4  K  I  It  4 
IF(K«I)  2*  4t  2 

2  C  •  A(Kt  I) 

DO  9  L  •  It  4 

9  A(Kt  L)  ■  A (Kt  L)  •  C  •  A< I t  Li 
AIK*  I)  •  •  C  /  B 
4  CONTINUE 
9  Allt  I)  ■  1.0  /  B 
RETURN 
END 
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•—n—  r. 


‘  -Ju*  y  w~1»l 


PAGE001 


AUTOMATH  SYSTEM  06/22/64 

•  •  *  SOURCE  PROGRAM  LISTING  *  *  * 


SUBROUT I NECOE 

COMMON  GNU*  RADO*  RADN,  NR*  NP«  PMU*  wO*  H*Xl *21 *X2 .22 . 

1S1*  52,  S3*  54,  F*  FF,  G,  GG,  SN«  ST 

DIMENSION  Xl<4) ,21 <4) ,22(4)  ,X2(4) ,  51(4*51  )«S2(4*5l  ), 53(4.51), 

1S4 (4,51  )•  F  (51) «  FF (51 ) »  G(S1).  GG (51),  SN(51)«  ST  <31 1 
1 »TRX (4,4) *  A (4) 

NRR  ■  (NR/NP) *1  1 

00  •  GNU/ 1 1.0  •  GNU)  2 

DO  10  I  •  1,4  3 

TRX(l.I)  ■  SI (1*1)  *00  •  (52(1 ,1) /RADO  ♦  0.3*54  ( 1 , 1 ) )  4 

TRX(2,I)  a  SKI, NRR)  ♦  DO  •  <S2(!«NRR)/RADN  ♦  0.5*  54(1, NRR))  5 

TRX (3* I )  a  (2, 0/7. 0)*S3(I.l)  -  52(1.1)  6 

TRX(4«!)  •  ( 2.0/7. 0) *53 ( I *NRR)  •  52(1, NRR)  7 

10  CONTINUE  8 

CALL  INV (TRX)  9 

00  20  K  al,4  10 

A(K)a-MO*(TRX(K«l)  ♦  TRX<K*2) )*1. 25*00  11 

20  CONTINUE  12 

So  30  IJ  a  1,  NRR  13 

F ( I J)  a  A( l)*S2 ( 1* I J) *A<2) #52 (2,1 J) *A (S) *S2 ( 3. 1 J) «A (4) *S2 (4,1 J)  14 

G( IJ)  a  A ( 1 )  *S4 ( 1 , 1 J) *A (2) *S4  (2, 1 J) *A (3) *54 (3  *  I J)*A (4) *54 (4, 1 J)  15 

FF(IJ)  *A(1)*S1 (1 . IJ) *A(2)*S1  (2*IJ)*A(3)*S1 <9* I J)*A(4)*Sl (4,1 J)  18 

GG(IJ)  *A(l)#SS(l.I J) *A (2)  *53 (2,1 J) *A(3)*S3 (3, 1 J) *A (4) *SS (4,1 J)  17 

30  CONTINUE  II 

RETURN  19 

END  20 
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AUTOMATH  SYSTEM 


04/22/44 


PAGEOOl 


•  ••  SOURCE  PROGRAM  LISTING 


SUBROUT I NESTRESS 

COMMON  GNU*  RAOO.  RADN.  NR*  NP*  PMU*  WO*  H.Xl.Zl .X2.Z2. 

1SI*  S2 •  S3.  34*  P*  PP.  G«  GG«  SN*  ST 

DIMENSION  XI  (4) *21 (4) *22(4) «X2 (4) *  31(4*31  1*12(4.31  ). 53(4, 51). 
1S4 (4.31  >.  P (31) .  PP(31) •  GtSl) «  GGI31) «  SN (51) «  ST (51) 

EE  «(1.0  •  GNU)/ (1*0  •  2*0*6NU) 

NRR  0 (NR/NP) *1 

DO  10  1J  •  1  .NRR 

SN(IJ)  0  2*0  *  PMU  •  EE  •  (WO  •  2*0  •  G(IJ>> 

ST(IJ)  0  -2*0  •  PMU  •  P(IJ) 

10  CONTINUE 
RETURN 
END 
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APPENDIX  D 


PURE  TORSION  IN  LARGE-DEFORMATION  ELASTICITY 
OF  AN  INCOMPRESSIBLE  ELASTOMER 


Take  as  a  reference  frame  cylindrical  coordinates  r,  0  ,  z  in  the  deformed 
body.  Then  a  point  r,  0  ,  z  was  initially  at  the  point  p»  9-f  z »  z .  where  p  is 
assumed  to  be  a  function  of  r  only,  and  the  point  of  the  undeformed  body 
is  given  by 


x-^  =  P  cos(9-Cz), 

x2  =  P  sin(e-Cz), 


x^  =  z, 


(D-l) 


where  £  is  the  angle  of  rotation  in  the  z  direction. 

The  components  Gjj,  G1J  of  the  metric  tensor  of  the  deformed  body  are 


f 
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The  components  g--,  g 

ij 
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metric  tensor 
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undeformed 
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0 

-Cr2 
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1 

(D-2) 

& 

(D-3) 


(D-4) 

&£ 

(D-5) 


2 

Since  G  =  g  =  r  ,  the  incompressibility  equation  =  1  is  satisfied. 


Consider  the  Mooney  stress-strain  relation  for  which  the  strain  energy 
functions  are  given  in  terms  of  two  constants,  Cj  and  Cz- 


Thus,  from  equations  (E-7)  and  (E-7a), 


i  =  2C, ,  Y  =  2C, 


«ik  = 


2  2 

2+Cr 


0 


0 


0 


l  0 


c2Hj  e 

r 

C  2 


(D-7) 


The  contravariant  stress  tensor  is 


11 

T  = 


2(ci  +  2C2)  +  2C2  C2  r2  +  P 


2  22 

r  t 


2  2 


=  2(C1  +  2C2)  +  2(C1  +  C2)  £  r 


+  P 


33  _ 


T  = 


2(ci  +  2C2)  +  P 


.33  = 


=  2(C1  +  C2)C 


x*- 


T12  =  0 


(D-  8) 


The  equilibrium  equations  are 


_ik  ,ri  rk,rk  ir  _  n 
t  .  +  r.  t  +  r.  t  =  0. 
l  lr  lr 


(D-9) 


From  the  metric  tensor  for  a  deformed  body, 


r  = 
22 

-2 


-r 

2 


r *  =  =  i/r ' 

1 12  21  ' 


(D-  10) 


The  rest  are  zero. 


89 


The  equilibrium  equations  are 


or 


From  equation  (D-7), 


|f  +  (4C2  -  2Cl)C2r  =  0 

and 

P  =  (^  -  2C2)C 2r2  +  a. 

Since 


(D-ll) 


(D-  12) 


(D  -  1 3) 


t^"  =0  at  r  =  r  (outer  radius) 

o 

P  =  -2(CX  +  2C2)  -  2C2C2ro2 

a  =  -2^  +  2C2)  -  ClC2ro2  (D-  14) 

P  =  -2(CX  +  2C2)  -  ci^r02  +  (Ci  “  2C2)C2r2  (D-15) 
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tU  =  -CjC2!^2  +  CxC2r2  =  C^r2  -  rQ2)  ^ 
r2r22  =  301C2r2  -  ^C2^2  =  C1£2(3r2  -  rQ2)  > 
t33  =  -C-,C2ro2  +  (Cx  -  2C2)£2r2 

t23  =  2(ni  +  C2)£  J 


The  physical  components  of  stress  are 

°ij  =  ^  (no  sum) 


C11 

II 

3 

k- 

ii 

c1C2(r2  -  ro2)  \ 

a22 

=  A22 

=  C1£2(3r2  -  ro2)  ^ 

ct33 

“  T33  = 

C1C2(r2  -  r02)  -  2C2C 2r2 

°23 

=tt23  = 

2(C1  +  C2)rC  / 

°31 

=  012  = 

0 

The  surface  tractions  at  plane  end  (bonded  surface)  are 


T  =  0 
r 


22 


Te  =  rT 

T  =t33 
z 


(D-  16) 


(D-17) 


(D-  18) 


(D-  19) 
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The  moment  required  for  the  twist  is 

M  =  J  (rTS)rdrde  =  J  r3T22drd0 

=  2n  [*  C1C2(3r3  -  r  2)dr 
J  b 


,  2  2 

4  r  r  r 


=  2nClC2  -  -  Vl,"  -  ^ClC2  t+  -  -t 

H  =  |  ClC2  [r04  -  r±4  +  2r02r.2]  . 


4  2  2 
r.  r  r. 

4  2  J 


(D-  20) 


Force  is  also  required  to  prevent  elongation: 


,b  b 

N  =  2 n  J  T33rdr  =  2tt  /  [C  £2(r2  -  r  2)  -  2C  £2r2]rdr 

a  -'n  °  2 


=  2n{-CiC2ro2/2  (rQ2  -  r.2)  +  -  2C2)C2  UrQU  -  r.4)} 

=  ^r2(  «  4  ,  fl  22  f*  ,  4  Us, 

^  l"  4  ro  T  ro  ri  2  (ro  ™  ri »  • 


(D-21) 


It  is  found  that  in  order  to  maintain  this  state  of  deformation,  pressure 
must  be  applied  at  the  inner  curve  surface.  The  magnitude  of  the  pres¬ 


sure  is 


-c1c2(ro2  - 1^2) 


where 


^  =  inner  radius 


0  =  outer  radius, 


APPENDIX  E 


LARGE-DEFORMATION  ELASTICITY  THEORY 
FOR  AXIAL  COMPRESSION 


In  order  to  characterize  the  deformation  of  an  elastic  body,  the  position 
of  each  point  in  the  body  before  and  after  deformation  must  be  described. 
That  is,  if  a  point  in  the  undeformed  body  is  at  y^,y^»y^»  then  in  the  de¬ 
formed  body  it  will  occupy  the  point  x^,x^,x^.  For  the  current  problem, 
it  is  convenient  to  employ  cylindrical  coordinates  for  this  description; 
therefore,  an  undeformed  system  (r ,  9'  ,z'  )  is  defined  by 


=  r  cos0 ’ , 
and  a  deformed  system 


1 

x  =  p 


COS0, 


y2  =  r  sin© ' , 
(p,9,z)  by 

x2  =  P  sin©, 


(E-l) 


(E-2) 


Tn  particular,  since  the  loading  in  the  present  problem  is  axially  sym¬ 
metric,  the  deformation  may  be  considered  in  the  form 

P  =  p  (r,  z 1 ) ,  z  =  z(r,z '  )  and  0  =  S'  +  tCr^O.  (E-3) 


In  order  to  determine  these  functions,  it  is  useful  to  consider  another 
description  of  the  deformation.  For  this  purpose,  the  distortion  which 
the  initial  polar  coordinate  system  undergoes  as  the  body  deforms  can  be 
taken  into  account  if  these  coordinates  are  considered  to  be  embedded  in 
the  material.  To  characterize  this  distortion,  the  metric  tensors  Gjj  and 
g-j,  which  are  associated  with  the  polar  coordinates  before  and  after  de¬ 
formation,  are  defined  as 


G  =  &L.  -&L- 

dcp1 


and 


g 


r  *  r 
Bx  ox 

B91  a 


(E-4) 


where  cp^  =  r,  =  CfP  =  z'.  Using  equation  (E-l) 


\xr  nVit  a  i  n 


(E-5) 
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and  from  equations  (E-2)  and  (E-3), 


gij 


r/SP\2  +  n2  /M,n2  .  ,dzx2-,  2  8t 

r(d?  +  p  (a7}  +(d7)]  p 


r^£.  M  +  02  at  iL_  +  L-  Sz_] 

^  d  r  d  Z '  Sr>  ^7*  1 


dr  dz'  dr  dz1 


\ 


p2  St 
p  dr 


rSP  dp  ,  2  dt  dt  ,  £z  dz_-i  2  dt 
Ldr  dz'  P  dr  +  ~J  P  ~ 


dr  dz  dr  dz’ 


dz'  LVBz 


2  dt_ 
dz' 


^^T)2  +  p2  + 


^z'J 


(E-6) 


In  terms  of  the  above  metrics,  Gjj  and  g^j,  the  stresses  t1^  in  the  body- 
may  now  be  determined.  As  a  first  approach  to  the  solution  of  the  prob¬ 
lem  under  consideration,  we  will  consider  the  material  to  be  incompres¬ 
sible  and  use  constitutive  relations  given  by 


-ij  = 


=  $Gaj  +  YQ1J  +  g 


ij 


(E-7)* 


whe  re 


AJ  =  r.rs 


g_.  G 
rs 


^  -  Gir  r,is 


g 


rs 


(E-7a)* 


Y  and  $  are  constants,  and  p  is  a  scalar  function  of  position.  Further¬ 
more,  from  equations  (E-5)  and  (E-6),  the  contravariant  components  of  the 
undeformed  and  deformed  coordinates  which  appear  in  equation  (E-7)  are 
found  to  be 


and 


% 


(E-8) 


Matrix  where  each  component 
is  the  cofactor  of  the  cor¬ 
responding  component  in  array 
.  (E-6).  This  is  shown  later  as 
\(E-21). 


\ 


/ 


(E-9) 


*See  reference  2,  equations  1.  16.  13  and  1.  15.  9,  page  28. 


Since  the  stresses  are  defined,  it  is  now  possible  t  obtain  a  set  of 
equations  whose  solution  will  yi  Id  the  desired  functions  (E-3).  These 
equations  are  obtained  by  usmg  the  equilibrium  equations 


T 

> 


id 


i 


Trj  +r>  j  Tim  =  0 

lm 


(E- 10) 


where  the  st  are  the  Christoffel  symbols  of  the  second  kind  defined  by 


rstr  “  *  Skr  [«sk,t  +  gtk,s  •  g5t,k] 


(E-ll) 


and  the  stresses  are  given  by  equations  (E-7). 


It  may  be  observed  that,  due  to  the  complexity  of  the  tensors  g*J  and  gjj, 
the  equilibrium  equations  as  defined  above  will  be  extremely  cut  ibersome 
(see  Appendix  D  equations).  It  may  also  be  observed  that  if  we  had,  in 
relation  (E-4),  identified  cp^s  with  the  deformed  rather  than  the  undeformed 
coordinates,  the  resulting  form  of  the  equilibrium  equations  would  have 
been  greatly  simplified.  The  following  considerations,  however,  indicate 
that  this  approach  is  not  practical,  since  it  introduces  difficulties  in  the 
boundary  conditions.  If  the  undeformed  coordinates  appear  as  the  inde¬ 
pendent  variables  in  the  equilibrium  equations,  as  is  currently  the  case, 
the  boundary  is  located  by  specifying  its  undeformed  position.  If,  however, 
the  equilibrium  equations  are  simplified,  as  described  above,  the  de¬ 
formed  coordinates  will  appear  as  the  independent  variable  and  the  boundary 
must  be  located  by  specifying  its  deformed  position,  which  is  unknown 
until  the  problem  has  been  solved.  It  is  observed  that  the  three  equations 
(E- 10)  contain  four  unknowns:  p,t,z,  and  P.  However,  using  the  fact 
that  the  material  is  incompressible,  a  fourth  equation  is  obtained, 


|G.  .| 


g 


ij 


=  1  , 


(E- 12) 


which  states  that  the  volume  of  any  portion  of  the  body  remains  constant. 


Up  to  this  point,  the  loading  history  of  the  body  has  not  been  discussed. 

However,  it  is  apparent  that  unlike  in  a  linear  elastic  theory,  the  order 

of  loading  must  be  specified.  In  order  to  build  this  information  into  the 

equations,  the  following  considerations  are  used:  Let  p,t,z,  and  P  be 

made  up  of  contributions  due  to  compression  and  torsion  (i.  e.  ,  let 

p  =  p  +  p  ,  z  =  z  +  z,  ,  etc.  ,  where  p  and  p.  represent  components 
C  t  C  t  C  w 
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after  compression  and  tension,  respectively).  Since  the  compression  is 
applied  first,  equations  (E- 10)  are  solved  for  pc ,  zc ,  and  pc(i.  e.  ,  tc  =  0  ) 
with  Pt  =  tc  =  zt  =  =  0.  Having  solved  for  the  compression  contribution, 

the  torsion  contribution  is  now  determined.  It  is  observed  that  the  torsion 
considerations  are  influenced  by  the  initial  compression  by  wa>  >f  the 
compression-torsion  coupling  terms  which  appear  in  the  equilibrium 
equations. 


Having  set  up  the  field  equations,  we  must  now  consider  the  boundary 
conditions. 


Figure  15.  Sketch  of  Coordinates  for  One  Lamination. 


For  the  compression  problem,  the  boundary  conditions  which  must  be 
satisfied  are  as  follows: 

at 

z'  -±h,  zc“+(h-6o),  Pc  =  r  (E-13) 


at  r  =  R0,  R^,  traction  is  zero,  or  in  terms  of  a  coordinate  system  on  the 
boundary,  Tnn  =  Tns  =  0, 


or 


(E-  14) 
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And  for  the  torsion  problem, 


at 


*c±(h-«o)>  pt=Pc 
(specified  twist  angle) 


r 


at 


r  =  R  ,  R. 


nn 


=  T 


nn 


0,  t. 


ns 


-  T 


ns 


=  0  and  t 


e 1  e  > 


=  o 


Having  formulated  the  problem,  a  method  for  obtaining  the  solution  will 
now  be  considered.  For  this  purpose,  consider  rewriting  the  equilibrium 
equations  and  boundary  conditions  using  the  displacements 


P 


c 


Zc  “  z1  wc>  “  z*  etc., 

rather  than  the  deformed  positions,  as  dependent  variables. 

If  it  is  considered  now  that  the  linear  elastic  equations  are  given  in  the 
form 


f(u  ,  **P)  =  0,  g(u,w,p)  =  0,  h(u,w)  =  0  , 


then,  with  the  above  change  of  variables,  the  nonlinear  equations  will  be 

f(u,w,p)  =  :;'(u,w,p),.  g(u,w,p)  =  G(u,w,P) 

and 


h(u,w)  =  H(u,w) 

where  F,  G,  and  H  represent  nonlinear  contributions.  The  solution  to 
these  equations  can  now  be  obtained  by  using  a  numerical  technique 
similar  to  that  used  in  the  solution  of  the  linear  equations.  To  do  this, 
the  nonlinear  terms,  F,  G,  H,  which  appear  in  the  equations  are  treated 
as  if  they  represent  a  nonhomogeneous  portion  of  the  linear  elastic  equa¬ 
tions.  That  is,  in  each  iteration  it  is  considered  that  the  nonlinear  terms 
are  known,  the  value  being  given  by  the  value  obtained  from  the  previous 
iteration. 
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Combined  Compression  and  Torsion 


Assume  that  a  rubber  annulus  is  placed  under  axisymmetric  compression 
followed  by  torsion.  Consider  deformation  described  by 


deformed  coordinates 
undeformed  coordinates 
where 


(P  ,0,z) 
(r,6',zr) 

P  =  p (r,z ' ) 
z  =  z(r,z’) 
0  =  0'  + 


The  metric  of  cp  in  the  undeformed  coordinates  is 


_  dyr 


where 


and 


1 

y 

2 

y 

T3 


r  cos©  1 
r  sinQ ' 


G.,  . 

lJ  Sep1 


cp  = 


V3 


r 

0T 


z'; 


then 


(E-  15) 


G.  .  = 
ij 


(E- 16) 


(E-  16a) 
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and 


1 


0 


0 


G1,3  =  I  0  l/r*  0 


0  0  1 


IG1*3!  =  l/r2. 


The  metric  of  cp  in  the  deformed  coordinates  is 


dxr  dxr 


SiJ  =  a?1  acpJ 


where 


x  = 


X3  = 


p  cosQ 
p  sinG 

z  . 


Expanding  the  separate  derivatives  in  equation  (E-18)  yields 


dx  _ 
dcp'*' 


dx  _  dx  dp  ,  dx  dG  _  d£_  .  Q  dt 

dr  dp  dr  dG  dr  dr  dr 


M. 

dcp3 


dx2  _  dx2  dP  ,  dx1  dG 
dr  dp  dr  dG  dr 


sinG  +  p  cosG 

dr  dr 


dxl  = 
dcp3- 


dx3  dx3  dz  __  dz 
dr  dz  dr  dr 


dx~*~  _  dx3-  _  dx3-  dG  _ 
d92  db  SQ  dGT 


p  sinG 


dx  _  dx  _ 


dcp  dG 


'  dG 


dx2  dG 

T 


=  P  COS 


G 


d6 


^  = 


dcp*  dG 


=  0 
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(E-  17 ) 


(E-  17a) 


(E-18) 


then 


&L.  a 

.  s 

_sd 

2e_ 

+ 

S9  = 

cosS 

ae-.p 

sinG 

St 

Sep3 

Sz' 

SP 

Sz ' 

ae 

Sz ' 

Sz ' 

Sz ' 

Sxf  = 

=  M2 

- 

+  S XL 

M_  = 

sin© 

5E-  +  p 

cosQ 

£t_ 

Sep3 

Sz ' 

’  Sp 

Sz' 

S9 

Sz ' 

32' 

Sz ' 

a 

£>d  = 

=  a£s 

= 

£z_ 

= 

Sep3 

Sz ' 

Sz 

Sz' 

Sz'  ’ 

'r/SP.d  +  n2  4.  (h£\2l  n2  at  r££-  +  n2  ^  £L-  +  M  M-V 

+p  W  +  P  L^I77+P  *r  *,«  +  Sr  ;w.tJ 


'Sr' 

,2 

Sr 


Sr  Sr  Sz' 


Sr  Sz 1 

,2  St_ 

Sz 


9L  +  p2M3L  +  iz  Sz_,  2  £t_  [(^-)2  +  p2  (£L-)2  +  (|£-)2] 

^  Sr  Sz1  H  >nl  ^  H  l's-ti'  k  v^^>/  ^7.'  J 


Sr  Sz’  Sr  Sz* 


Sz1  Sz 


Sz* 


Sz1  , 

(E- 19) 


which  implies 


IgJ  =-P2  tr  ip2tr  [z2,+  p2t2,+  p2,]  -p\,  [prp7,+  p\t„+  z^,]) 


r  z 


r  z '  r  z 1 


+  p2  {[p2  +  p2t2  +  z2][z2  +  pVr+  p2  ]-[p  p  .  +  p2t  t  +  Z  Z  ][p  p 

I  X  XX  4  2  l  L  I  a  X  Li  X  £ 


+  p\t  ,  +  Z  z  ]}  -  p2t  .  {p2 t  ,[p2  +  P2t2  +  z2] 


r  z 


r  z ' 


z '  r 


-  P2trrprPz,  +  p  trtz,  +  zrzz  i  ] } 


=  P2[P2Z2,  +  Z2Z2,  +  P2P2,  +  Z^2,  -  P2P2,-  Z2Z2,  -  2PrPz,ZrZz,] 


r  z 


r  z 


and 


g 


ij 


=  p2  r  1^;-  -  =  rph(r,Z')]: 

Sr  Sz'  Sr  Sz' J 


(E-20) 
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Thus 


f2  rM  <£_  +  £2  +  it  rr^i-)2  + 

dz'  Ldr  dz'  dr  dz'  b r  dz1'  'bt' 


[Ph]2 


.2  ,^t_  rSfi.  ££_  4  la  *i_l  *  M  r/£i_)  +  f^ie _ )2ll  n2  I  (&■  -  &  &—)2  +  (&■  £t-  -  —  £i)21  +  (&■  £*-  . 

P  ldz'  ldr  3z'  dr  dz' J  drUdz''  lpz';j)  P  Udr  dz'  drdz';  ldr  dz'  dz'  dr'  J  'dr  dz' 


2  rde.  d£_  dj,, 
IT  dz'  dr  dz'J 


\ 


p2  r££_  r(££.\2  +  (£i)2]  +  &  fit  &-  +  if  £f_-]i 
P  ldz'  udr'  ldr;  1  dr  ldr  dz'  dr  dz,JJ 


.2  rjfi.  ££_  +  is  £l_-| 
dr  dz'  dr  dz' 1 


\ 


2  rdt  r/£B.\2  4  f£f)2i  4  £i  riB.  ifi—  4  if  £f-]l 
P  dz '  'dr'  dr  J  dr  'dr  dz'  dr  dz,J) 


P2  tf/  ♦  <j*M 


Now  the  strain  invariants  are 


I1  G  ^  gij  gll  +  ^2  g22  +  g33 


(E 


and 


|G.  .|  2 

t  -  i.r  _  r 

|  g.  .[  2,  2  * 

1 bij  P  h 


(E 


but  for  incompressibility,  1 3  =  1;  hence, 


P  = 


h(r,z’) * 


(E 
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»^sKw?r 


it 

dr  dz ' 


-21) 


-22) 


-23) 


-24) 


a 


Also,  from  equation  (E-7a), 


QiJ  =  Ix  G1J  -  Glr  GJS  grs 

( 1  °  °\ 

where  G1^  is  the  diagonal!  0  l/r  0 

\0  0  1 


Thus,  the  components  of  G1-1  are 


Thus , 


nll  _  T  rll  P11  rll  =  ,  /  2  , 

Q  Iq  G  —  G  G  §99  S 


^22  _  f  022  n22  n22  _  /  2 


>22  fa3  3 


I,  G^  -  G^  g22  =  l/r  (gxl  +  g33) 
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_  T  r33  r33  n33  =  _  +  l/r2 

I1  G  “  u  G  g33  gll  g22 


Q12  =  Qa  =  G11  G22  g12  =  -  l/r2  gi2 


Q13  =  Q31  =  -  G11  G33  gy  =  -  gy 
q23  =  q32  =  .  g22  G33  g23  -  -  l/r2  g 


>23 


t1j  =  $  G1'-’  +  Y  +  P  ^ 


deformed  undeformed 
coordinates  coordinates 


defornled 

coordinates 


/•  (E- 25) 


J 


) 


which  is  then  substituted  into  the  equilibrium  equations,  (E- 10). 


Compression  Only 

For  the  problem  of  axial  compression  only,  the  formulation  is  carried 
further  as  follows: 
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For  the  case  of  notation,  let 


From  equation  (E-17), 


0 

1/r2 

0 


From  equation  (E-20), 


h  =  fM  dz  _  _5z  Bp  > 
or  BzF  Br  3zF^ 


[u  +  |J)a  + 


3z'; 


cw  Bu  i 
Br  3z,J: 


and  from  equation  (E-23),  h  is  also  given  by 

h  =  r/p  . 


From  equation  (E-21),  the  components  of  g1.)  are 


1!  =  JL  r(^-)2  +  (^-)2]  =  —  r(i  +  ^-)2  +  (^-)2^ 

g  .2  LVBz,;  ^Bz,;  J  .2  IU  Bz,; 
h  h 


‘3z 


22  &  1  rB£  £z_  £z  3£_-,2  =  1_  =  -1  =  hf 

9  9  1 >-n  >f7  I  “  ** 1  -1  2  2  2 

P  r2(1+U) 

r 


g  p2  ^ 2  1  Br  BzF  Br  Bz 


*33  .  i  r(M)2  +  (St/-,  =  -L  r(1  +  iu)2  +  (^[)2] 

g  u2  U3r'  ''Br'  J  .2  1  u  Br'  VBr' 

h  n 


13  =  -31  =  1_  r2£  S£.  Bz  Bz__1=  lr/,  ,  %2L.  +  &  (i  +  2k 

g  g  .2  L3r  Bz1  BrSz,J  ,  2LU  Br'  Bz»  Sr  u  Bz 

n  h 


12  _  21  _  23  _  _32  _  n 

g  =  g  -g  -  t  -  0 


(Since  torsion  is  zero,  •’  :  0) 
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'  •’  •■*■’** 


(E-  26) 


(E  -  27 ) 


(E-27  a) 

\ 


>• 
r)] 


(E-28) 


F'om  equation  (E-19),  the  components  of  g  —  are 


■v.  2  >2  >2  ^2 

=  (°L)  +  (M)  =  (X  +  221)  +  (-2H)  =  h2  g33 

gll  W  '3r'  u  3r;  v3r'  g 


l2  __  2 , 


g22  =  -  (u  +  r)2  =  r2(l  +  u/r)2  =  r2/h2 

g„  =  (|*7)2  +  (f^T)2  -  (1  +  (&)  =  h2  S11 


=33  vBz' 


oz 


3z 


oz 1 


=13  or  3z’  3r  dz 


Sr'oz'  dz 


dz 1 


£  ^  £  =S£3L  +  M1l=(1  +  £u)£u_  +  ^_(X  +  lw.)=  _h2Jl=h2B 

g31  g13  or  3z»  3r  dz'  U  dr'dz'  3z,U  dz'J  r  g 


s23  g32  g21  g12  0 


(E  -  29) 


Now  the  components  of  t1J  from  equation  (E-7),  by  substituting  from 
equations  (E-17),  (E-25),  (E-28),  and  (E-29),  are 


n  = ,  +  i_  +  [(i  + 1</  +  (m.)2]  (y  +  jl, 

h  n 


T22  =  l  + 1_  (  +  )  +  P  ic 

1  2  2  '*11  g33'  r  2 

r  r  r 


dU, 

dr' 


M 

dz 


rcW 

‘dr 


,  2 


=  $  +  Y  (g-n  +  l/r^  g22)  +  P  =  $  +  2 

h  h 


2  _  x  ,  „  gll  •  v 


=11 


+  rd  +  ^)2  +  (f*n  (y  +  p/h2) 


t13  .  t31  =  (y  +  £5)  [(1  +  f )  (^)  +  |  (1  +  &) ] 


dz 1 


t23  .  t32  .  t12  .  t21  .  0 


J 


(E-  30) 
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The  equilibrium  equations  are  now 


or 


-1JA  =  o 


Ti0  +  ri_  Trj  +  rj  in  =  0 
,  i  n  lm 


(E-H) 


as  shown  in  Equation  (E -  10), 
whe  re 

-r 


r 


st 


i  [Ssk,t  +  gtk,s 


"  g 


st,kJ 


(E-  32) 


Note: 


rr.  = 

st  ts 


In  the  current  problem  of  compression  only,  the  only  nonzero  Christoffel 
symbols  are 


ril  5  ^g  gn}i  +g  ^2gi3,i  “  gn,3^ 

,1 


r22  =  "i  (gl1  g22,l  +  g22,3} 


13 


r1  =  r1  =4  fg11  a 

13  31  5  ^g  gn,3 


+  g  3  633^) 


~\  Q 

r33  -  i  Tg  g33^  +  g  (Sg^^  -  g33jl) 

r2  r2  i  22 

12  r21  ^g  g22jl 


11 


p2  _  ^2  _  i  22 

23  32  g22^3 


rii  ^g  3  gll,l  +  g33  ^2g13,l  "  gll,3^ 
f22  =  ^  g22,l  +  ^  g22,3^ 


f33  ^  ^  g^  g33,3  +  ^  ^2g13,3  "  g33,l^ 

r31  =  ri3  =  ^  ^  3  sll,3  +  g33  g33,l"* 


31 
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Now  the  equilibrium  equations  become 

£  *  3  *  ’  rn’“  *  3  fi'3*  *  *  1/* 

*  r»,3‘  *  ru’u  *  ’ll  *  if! ' 0 


(E-  34a) 


and 


+  2  r33T33  +  3  r33T13  f  r32x33  +  riy13 
+  r^r13  +  r^x33  +  r3^11  +  r3^22  =  o  . 


(E-  34b) 


{Note:  The  quantities  underlined  are  those  that  contribute  linear  terms.  ) 

The  third  equilibrium  is  satisfied  identically,  as  expected,  since  there  is 
no  dependence  upon  angle  Q  (axisymmetric).  Now  consider  rewriting  the 
stress  component  in  the  form  of  linear  terms  plus  the  nonlinear  contri¬ 
bution.  To  do  this,  consider  first  the  continuity  equation  from  equations 
(E-27)  and  (E-27a). 


For  convenience,  the  following  notation  will  be  used: 


2*  = 


u 


2u 


2r  2z '  V*  etc* 

From  equations  (E-27)  and  (E-27a) 


(E-  35) 


(1  +  u  )  (1  +  w  ,)  -  w  u  , 
r  z  r  z 


or 


1  + 


(1  +  ur)  (1  +  wz)  (1  +  “  wr  2I 


(1  +  U)  =  1. 


Expanding,  this  yields 


(E-  36) 


u  +  w  +  —  =  —  (u  w  .  +  ^  w  ,  +  ^  u  ■+■  u  w,~  +  w  u  (l  +  ”) 
r  z1  r  r  z'  r  z 1  r  r  r  z1  r  r  z'  r 7 


or 


+  W  .  +  ^  =  A 


U„  T  w_, 
r  z 


(E-  37) 


where  A  =  -  ur  wz,  +  ^  wz'  +  r  ur  +  ur  wz!  r  +  wr  V  ^  +  r^* 


Let  (r,z)  be  initial  coordinates  and 
(r',z')  be  final  coordinates. 
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Equilibrium  equations  are  written  with  r  and  z  as  independent  variables. 

A  location  is  then  described  by  its  initial  coordinates;  one  then  solves  for 
r'  and  z'.  The  stresses  are  expressed  in  the  final  coordinates;  and  on  the 
traction-free  boundary,  the  free  inner  and  outer  edge,  the  normal  and 
shear  stresses,  o  nn  and  ans,  are  required  to  be  zero.  To  express  this 
in  terms  of  oriri,  a  z?z  1  >  and  C  z  1  r  1 ,  which  are  in  turn  functions  of  the  co- 


2k. 


ordinates,  one  needs  the  boundary  shape.  That  is,  the  quantity  QZ  is 
needed  where  the  boundary  location  is  given  as  r  =  R.  If  a  finite -diffe rence 


P  .  +,  -  P  '  • , 

method  is  to  be  employed,  the  derivative,  S£  ,  is  equivalent  to  - - : 

dz  z  -  z'  .. 

V  -L-  L 

With  this  approach,  the  grid  shape  is  constant  and  the  equilibrium  equa¬ 
tions  are  written  in  the  deformed  coordinates,  which  is  the  more  con¬ 
venient  form. 


Now,  if  we  let  P  =  P'  -($+24),  (E-  38) 

then  expressions  (E-30)  yield 

t11  =  {  +  Y  (1  +  *)  +  (1  +  2„z,  +  vz,“  +  uz /)  [V  +  (pt  -  {  -  2:)(1^)  J 

2 

t11  =  ($  +  y  +  2Y  £  +  ©  Y)  +  [Y  +  P'-$-2Y+2  £  P'-2  7  ($+2Y) 

2 

+  (P'-$-2Y)  (~) ]  +  2w  ,Y  +  2w  f  (u'-4-2: )+  2w  ,  (P'-*-2Y)  [2  7  +  ©  j 

Z  Z*  «  A  1 

2 

+  (w  ,2  +  u  ,2)  [Y  +  (P'-$-2Y)  (1  +  7)  ]• 
z  z  * 

T11  =  2Y  U  -  2  ^  ($  +  2Y)  +  2w  . ( $+2Y )  +  P'  +  B  (E- 39) 

r  r  z 


where 

2  2 

B  =  ©  Y+  2  £  P'  +  (P»-$-2Y)©  +  2wz,  (PT-$-2Y) [2  7  +  (~)2] 

2 

+  (w  ,2  +  u  ,2)  [Y  +  (p'-$-2Y ) (l  +  7)  • 
z  z  * 

Then 

t11  =  -2($  +  Y)(wz,  +  J)  +  P'  +  B, 
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and  from  continuity,  equation  (E-37), 

*  -K,  +  S)  =  -  A 


or 


also, 


11 


=  2($  +  Y)  u  +  P'  -  2($  +  Y)A  +  B 
r 

l - 1 — ; 

nonlinear 


(E-4 1 ) 


linear  terms  can 
identify  Lame 1  s 
constant 


contribution 


t13  =  -  [Y  +  (P'-*-2*)(l  +  7)  [uz,  +  wr  +  ur  v  +  wr  wz,] 

2  u  2 

=  _  {-($  +  y)(u  f+  w  )  +  [P'(l  +  £)  -  (*  +  2Y)(2  £  +  (“)  )] 


[uz'+wr  +  ur  uz'+ Wr  Wz,]  +  Y  [Ur  V  +  wr  wz']] 


(4  +  Y)  (u  ,  +  w  )  +  C 

-  j8  "  ^ 

linear  terms  with  nonlinear 

p,  =  $  +  Y  contribution 


(E-42) 


where 


C  = -l[P'(l+“)2  -  (4+2Y)(2  “  +  (;)  )]  Cu7.i+Wr+Ur  V+Wr  z'] 


+  YfUj.  uz,  +wrwa,]]j 


also, 


22 

T  = 


{$  (1+7)  +  y(i  +  7)  [z  +  z  u  +  a/  ,  +  U  +W  . 


r2  (1  +  z) 


+  w2  +  u2,]  +  P'-4-2Y) 
r  z 


r2  (1  +  Z) 


-i - -  [2  “  (4+2Y)  +  2Y(ur+wz,)  +  P'  +  2Y  <*) 


+  Y[2  ^  +  (^2][2V  ^  z-  +  wr  +  V]] 
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from  continuity  u  +  w  ,  =  A  -  or 
7  r  z '  r* 


r22  = 


-2  (1  + 


2  [2($  ty)  J  +  P'  +  2YA  +  D] 


(E-43) 


71 

linear 


nonlinear 


where 


2  2 

D  =  2Y  C11)  +  Y[2  -  +  (-)  ][2u  +  2w  ,  +  u2  +  w2.  +  w2  +  u2,]  . 

r  r  r  J  r  z 1  r  z 1  r  z ' 


It  is  seen  here  that  the  factor 


1 


r2  (1  +  “  )2 
r 


2  2 

causes  the  linear  terms  in  t 


22 


to  be  different  from  the  expression  for  t  encountered  in  a  linear  analysis, 
This  is  due  to  the  fact  that  the  physical  and  tensoral  (  T  ■' ^ )  components  of 
the  stress  tensor  are  different.  Their  contribution  in  the  equilibrium 
equations,  however,  will  be  the  same. 


Also,  from  the  similarity  in  the  form  of  and  t  ^ ,  the  following  is 

immediate]  /  obtained: 


,33 


+  Y)  w  ,  +  F1  -  2($  +  Y)  A  +  E 


(E-44) 


where 

2  2  2 

E  =  (-)  +  2  ~  P'  +  (P'-#-2Y)0  +  2u  (P'-*-2Y)[2  “  +  0  1 

'  J'  i  I  ^  A 

2 

+  (u  2  +  w  2)[Y  +  (P'-$-2Y)(l  +  ] 

'  r  r  r 

and  the  first  of  the  equilibrium  equations,  (E-34a  and  b),  yields 


2(M)  ^  -  2(i+f)  M  +  |S  +  (f+Y)(93i_  _£*_•>  +  3C_ 

3r2  Br  3r  dr  '  3zi2  3z' 


+  2  (♦*,)  *  +  £+[§»-  “][— ^][2U«)  *  +  p.] 


(1  +“)' 


dr 


■F  [1  +  f^lD - ”  “][-2  ($+y)A  +  B]  -l/r  [2($+Y)il  +  P>] 

r  (l  4  -)  r 

r 
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-  h  f[l  +  +  d  +  V  +  V  +  [2  &  +  <t^ 

r 

+  (§“,)  ][(l  +  “)(1  +  §£)]}  t2(tn')  “  +  P>  +  2T  A  +  D) 

-  Z.  [?;'  A+D]  +  — - i— r  g13  g  [2(i^Y)  “  +  P'  +  2Y  A  +  D] 

r  r  (1  +  r 


+  [other  nonlinear  terms  from  expression  E-34a  ]• 


This  equilibrium  equation  may  be  written 
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and  the  second  of  the  equilibrium  equations,  (E-37b),  yields 
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